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, The design and construction of a number of new and more sensitive detectors of grav- 

' itational radiation is currently underway. These include the LIGO detector being built 

in the United States by a joint Caltech/MIT collaboration [(Abramovici et al. 1992)| , the 
■ VIRGO detector being built near Pisa by an Italian/French collaboration (Caron et al. 

1995), the GEO-600 detector being built in Hannover by an Anglo/German coUabo- 
ration (Danzma nn et al. 1995), and the TAMA-300 detector being built near Tokyo 
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, (Tsubono 1995) . There are also several resonant bar detectors currently in operation, 

' and several more refined bar and interferometric detectors presently in the planning and 

, proposal stages. 

It's not clear when the first sources will be detected - this may happen soon after the 
detectors go "on-line" or it may require decades of further work to increase the sensitivity 
^ i' of the instruments. But it's clear to me that eventually, when the sensitivity passes above 

some threshold value, the gravity- wave detectors will find sources. 

The operation of these detectors will have a major impact on the field of gravitational 
physics. For the first time, there will be a significant amount of experimental data 
to be analyzed, and the "ivory tower" relativists will be forced to interact with a broad 
, range of experimenters and data analysts to extract the interesting physics from the data 

}^ ' stream. Being an optimist, I think that this "interesting physics" will include sources 

, of gravitational radiation which we have not hitherto expected or even conceived of. It 

promises to be an exciting time. 

There are at least three major types of "known" sources (Schutz 1991; Thorne 1995). 
These are (1) coalescing binary systems, composed of neutron stars and/or black holes, 
(2) pulsars (and other periodic sources), and (3) supernovae (and other transient or burst 
sources). These are listed, roughly speaking, in increasing order of "detection difficulty". 
For example, the waveforms (chirps) of the coalescing binary systems are already known 
to a high enough degree of precision (Blanchet et al. 1995; Apostolatos et al. 1994; 
Tagoshi & Sasaki 1994) to enable reliable detection, using the classical signal analysis 
technique of matched filtering (Dhurandhar & Sathyaprakash 1994; Apostolatos 1995). 
Some order-of-magnitude studies have shown that the amount of required processing can 



be done with a dozen or so high-performance desk-top workstations (Owen 1995)| . The 



major unknown for this type of source is their number density, which is uncertain by 
more than an order of magnitude. Hence the rate of events which occur close enough 
to us to be observable can not be estimated very precisely. The second type of source 
is pulsars. Once again, the form of the signal is known very precisely: it's just a sine 
wave in the solar system barycenter, a coordinate system at rest with respect to the 
Sun KNicbauer et al. 1993)| . In this case however, the signal processing problem is far 
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more difficult, because the orbital motion of the earth around the sun and the rotational 
motion of the earth around its axis modulates the pulsar frequency. This means one must 
separately analyze the signal for each of ~ 10^'' separate patches on the celestial sphere, 
each of which would have a distinctive pattern of frequency modulation. The required 
processing speeds are currently beyond the limits of even the most powerful computers. 
Of course, it may be just a question of waiting until a better search algorithm is developed, 
or faster computers become available! The third source which I've listed, supernovae, is 
perhaps the most difficult to detect. There are two reasons. First, we do not have precise 
predictions of the gravitational waveform such an event would produce and it's hard 
to design a data analysis algorithm to look for something unknown. Furthermore, these 
explosions are expected to be rare events and because certain kinds of instrument noise 
might look quite similar to supernovae, I chose to list them last, and categorize them as 
the most difficult sources to detect. 

The subject of these lectures is a fourth type of source, quite different in character from 
the three listed above. These arc the "stochastic" or "background" sources (Michelson 
1987; Christensen 1992; Flanagan 1993). Roughly speaking, these are "random" sources, 
typically arising from an extremely largo number of "unresolved" independent and un- 
correlated events. This type of background could be the result of processes that take 
place very shortly after the big bang, but since we know very little about the state of the 
universe at that time, it's impossible to say with any certainty. Such a background might 
also arise from processes that take place fairly recently (say within the past several billion 
years) and this more recent contribution might overwhelm the parts of the background 
which contain information about the state of the early universe. 

These sources are "unresolved" in the following sense. If we study an optical source, 
somewhere in the sky, using a telescope with a certain angular resolution, then details 
of the source can be "resolved" if the angular resolution of the telescope is smaller 
than the angular size of the features or objects being studied. In the case of the LIGO 
experiment, and similar detectors, the angular size of the antenna pattern is of order 90°. 
Hence almost any source is "unresolved" in that it makes a significant contribution to the 
detector output for almost any orientation of the detector and the source. When many 
such sources are present, even if they are pointlike, the resulting signal has a stochastic 
nature. 

My motivation for studying these stochastic sources is two-fold. The first reason is 
a rather hopeful one. Because the gravitational force is the weakest of the four known 
forces, the small-scale perturbations of the gravitational field decouple from the evolution 
of the rest of the universe; at very early times. Currently, our most detailed view of the 
early universe comes from the microwave background radiation, which decoupled from 
matter about 10"'' years after the big bang, and gives us an accurate picture of the universe 
at this rather early time. Some rather simple estimates (which I'll elaborate later) show 
that if the current crop of gravity- wave detectors do detect a background of cosmological 
origin, then it will carry with it a picture of the universe as it was about 10~^^ seconds 
after the big bang. This would represent a tremendous step forward in our knowledge, 
and is the main reason for my interest. 

The second reason for my interest in stochastic sources is rather more practical. One 
often hears it said that searching for signals in the output of a gravitational wave antenna 
is like searching for a needle in a haystack. And indeed, for the first three types of source 
listed above, this is true. One has to search through tons of rock (the data stream) in 
order to find the one precious gem (say, a binary chirp) . Most of the rock is barren, and 
the challenge is to isolate the one tiny volume containing the material of interest. As 
I have already discussed, for gravitational wave sources, the amount of computational 
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power required for this careful search can be very large. However the situation is rather 
different when one is searching for a stochastic background. The analogy in this case is 
mining aluminum, where on average every ton of ore contains a certain number of kilo- 
grams of aluminum. The situation is analogous for stochastic gravity-wave sources. As 
you will see, whenever two detectors are operating simultaneously, even if only for a few 
seconds, we get a little bit more data and information about the stochastic background. 
And as you will also see, it is easy to analyze this data. The "signal" in this case is a 
very low bandwidth one, so the essential part of the data analysis for a stochastic back- 
ground can be done on a garden-variety personal computer. For example, in the case of 
the LIGO detectors, the part of the signal carrying a significant amount of information 
lies below a few hundred Hz (see Section 3^ and Fig. for details). Thus the rate at 
which information needs to be processed is only a few hundred data points/second; a 
very manageable rate. 

These lectures are organized as follows. In Section |^, I discuss some of the general 
properties that a stochastic background of gravitational radiation might have. I show 
how such radiation is characterized by a spectral function, discuss some of its statistical 
properties, and show during which cosmological epoch the radiation which falls into 
the bandwidth of the ground- and space-based detectors was produced. In Section || 
I show how one can combine data from two or more gravity-wave detectors to either 
put limits on the amplitude of a stochastic background, or to actually detect it. I 
do this in several steps, first giving a crude argument which demonstrates the main 
detection strategy, then discussing the reduction in sensitivity which comes about from 
the siting of the detectors and their relative orientations on the earth. This is followed 
by a rigorous derivation of the optimal signal processing strategy, and a calculation 
of the expected signal-to-noise ratio and the minimum detectable energy-density in a 
stochastic background. In Section |^ I discuss the observational facts: what we actually 
know about the stochastic background. There are strong limits on the spectrum of a 
stochastic background coming at very long wavelengths from observations of the isotropy 
of the Cosmic Microwave Background Radiation, and limits at higher frequencies from 
observations of timing residuals in millisecond pulsars. Finally, there is a limit on the 
integrated spectrum arising from the standard model of big bang nucleosynthesis. The 
last part of these lectures is far more speculative. In Section ^ I discuss some of the 
potential ways in which a stochastic background might arise from processes that take 
place early in the history of the universe. I discuss in some detail three particular models 
(inflation, cosmic strings, and bubble formation in a first-order phase transition) each 
of which gives rise to a different spectrum of radiation, and examine in some detail the 
mechanisms at work in each case. In each case I indicate if and when the experiments 
being built might detect these potential sources. This is followed by a short conclusion, 
and an appendix containing a few calculational details. 



2. The stochastic background: spectrum & properties 

2.1. Digression - the cosmic microwave background radiation 

In order to discuss the spectrum of gravitational background radiation, we need to in- 
troduce notation. I will do this with an analogy and an example. 

Let's begin by considering the electromagnetic background radiation, conventionally re- 



ferred to as the Cosmic Microwave Background Radiation (CMBR) (Kolb & Turner 1990) 



This radiation was originally produced when the universe had a temperature of about 
3000 K, at a redshift of approximately Z = 1100. Today, this electromagnetic radiation 
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has a Planck blackbody spectrum with a temperature of about T = 2.73 K. Each cubic 
centimeter around us contains « 400 CMBR photons of energy « 10~^^ erg. The total 
energy density in this radiation field is pem = 4.2 x 10"^"^ ergs/cm"^, and the characteristic 
frequency of the photons that make it up is (a few times) /o = kT/h = 5.7 x 10^° Hz. 
Here /i = 6.6 x 10"^'' erg — sec is Planck's constant, and fc = 1.4 x 10~^^ erg/Kelvin is 
Boltzmann's constant. 

In order to characterize the spectral properties of this radiation, we'll now consider 
how this energy is distributed in frequency. In a spatial volume V, the amount of energy 
dE contained in this radiation field between frequencies / and f + df can be expressed 
as 

..MW(;^)(^^). (.1) 

where c = 3 X 10^° cm/scc is the speed of light. The different factors appearing on 
the right hand side of this equation are (1) the number of polarizations (2) the energy 
per quanta (3) the number of quanta per mode, and (4) the number of modes in the 
frequency interval. Dividing both sides of this equation by the volume V, one may write 
the energy density within the frequency range df as 

^dE _ Snh fdf 

dpem = — - ^^^/^j, _ ^ . (2.2) 

It will prove convenient to write the differential energy density within a unit logarithmic 
frequency interval. This is 

dp.^ ^fdp^^ 3.8 ^ io-i4£IgJ (iVf _1Z1\ . (2.3) 



dlnf ' df ■ cm3 V/o/ 

This formula contains complete information about the spectral distribution of energy 

in the CMBR, and if wc were interested in discussing electromagnetic backgrounds, we 
would probably stop here. However there is a slightly different standard convention used 
to describe gravitational wave backgrounds. 

In describing gravitational wave stochastic backgrounds, it is conventional to compare 
the energy density to the critical energy density Pcriticai required (today) to close the 
universe. This critical energy density is determined by the rate at which the universe is 
expanding today. Let us denote the Hubble expansion rate today by 

Ho = hwo 100 = 3.2 X IQ-^^hwo— = 1-1 x IQ-^^chwo — - (2.4) 

sec — Mpc sec cm 

Because we don't know an accurate value for Hq (a matter of considerable controversy 

in the literature) we include a dimensionless factor of /iioo which almost certainly lies 
within the range 1/2 < hioo < 1- The critical energy-density required to just close the 
universe is then given by 

3c2ij2 

^critical = -^-^ « 1.6 X lO-'^hfoo crgs/cm^. (2.5) 

This leads to our fundamental definition in this section, of a quantity which we will be 
using for the remainder of these lectures. 

In the remainder of these lectures, we will be discussing gravitational wave, rather 
than electromagnetic backgrounds. However, to complete this section, we first define the 
analogous quantity for electromagnetic backgrounds. This quantity is a dimensionless 
function of frequency 

^ / ,-x 1 dpcm „x 




Frequency f (Hz) 



Figure 1. The solid curve is the fractional energy density r2em(/) contained in the 2.73 K 
electromagnetic background radiation (with hioo set to unity). The dashed curve shows the 
corresponding quantity for a 0.9 K blackbody. If the gravitational perturbations had been 
in equilibrium with the matter fields, this is the expected spectrum of the gravitational wave 
stochastic background. Both axes are logio. 



For the thermal spectrum of electromagnetic radiation in the CMBR, we have 



Jo J \ef/f«-l, 

A graph this function is shown in Fig. |l| Also shown in Fig. |^ is the spectrum that 
the gravitational wave stochastic background would have, if at early times in the history 
of the universe the fluctuations in the gravitational field had been in equilibrium with 
the other matter and radiation in the universe. In this case, the gravitational wave 
stochastic background would have a thermal spectrum, with a temperature of about 
0.9 K (Kolb & Turner 1990), It is smaller than the temperature of the CMBR because 
in a conventional hot big bang model, the gravitons would have decoupled when the 
temperature of the universe dropped below the Planck temperature, when the number 
of entropy degrees of freedom was 106.75 in the standard GUT model. Since the number 
of degrees of freedom today is only 3.91, the graviton temperature is less than that 
of the CMBR by the ratio (106.75/3.91)^/^ (See discussion in Kolb & Turner (1990), 
pg 75, between (3.89) and (3.90)). However it is unlikely that this equilibrium could 
have been established; the time required to establish the equilibrium is longer than the 
characteristic expansion time (the Hubble time) of the universe because the gravitational 
interaction is so weak. While it is therefore unlikely that this 0.9 K thermal spectrum is 
present, it is nevertheless a useful benchmark for comparison. 

A number of graphs similar to Fig. ^ will appear in these lectures, so a couple of 
comments are in order. First, the reader will notice that the horizontal axis encompasses 
an enormous range of frequencies. The lowest frequencies / « Hq are those of waves 
which only oscillate a single time in the entire history of the universe, and whose period 
is a Hubble time Hq^\ The highest frequencies shown are those of visible light. Second, 
these graphs make it easy to see the total amount of energy contributed by the radiation 
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to the energy density of the universe. From the vertical axis, one can immediately see 
from the graph of r2cm(/) that the CMBR contains, in the vicinity of 10^"'^ Hz, about 
10"^ of the energy required to close the universe. 

2.2. Notation - the spectral function ri(/) for gravitational waves 

In order to characterize the spectrum of a stochastic gravitational wave signal, we in- 
troduce a quantity for the graviton background which is analogous to rJom- The precise 
definition is 

n{f)^n,Af)^^—^ (2.8) 

Pcritical " in / 

The subscripts "gw" for "gravitational wave" are omitted when there is no danger of 
ambiguity. 

There appears to be some confusion about in the literature. Some authors assume 
that is independent of frequency /; this is true for some cosmological models, but 

not for all of them. The important thing is that any spectrum of gravitational radiation 
can be described by an appropriate ^{f)- With the correct dependence on frequency / it 
can describe a flat spectrum, or a black-body spectrum, or any other specific distribution 
of energy with frequency. 

You will also notice that it follows directly from the definition that the quantity 
^ioo^(/) is independent of the actual Hubble expansion rate. For this reason, we will 
often focus attention on that quantity, rather than on alone. 

It is sometimes convenient to discuss the spectrum of gravitational waves and the 
sensitivity of detectors in terms of a characteristic "chirp" amplitude. This is the dimen- 
sionless gravitational- wave strain h = AL/L that would be produced in the arms of a 
detector, in a bandwidth equal to the observation frequency. This is related to ft{f) by 

-100 Hz 



h,{f) = 3 X io-2"/iioovW)— (2.9) 
(Equation 65 of Thorne (1987)). Apart from an overall factor, this formula can be easily 



derived by dropping the time derivatives from (A 13). For example if 0(/) = 10 over 
a bandwidth 50 Hz < / < 150 Hz then the strain in an ideal detector is he ~ 10"^'*. 

2.3. Assumptions about the stochastic gravitational-wave background 

Does il{f) contain all information about the stochastic background? The answer is "yes" 
provided that we make enough additional assumptions. 

We will assume from here on that the stochastic gravity wave background is isotropic, 
stationary, and Gaussian; under these conditions it is completely specified by its spectrum 
ri(/). Each of these three properties might or might not hold; before moving on, let's 
consider each in turn. 



It is now well established that the CMBR is highly isotropic (Kolb & Turner 1990). In 
fact this isotropy is surprising, because in the standard model of cosmology, the angular 
size of the horizon at Z ^ 1100 is only about 2°. Nevertheless, it is experimentally 
well-established that the largest deviation from the isotropy of the CMBR arises from 
our proper motion with respect to the rest frame of the universe, at the level of 1 part in 
10^. The next largest deviations from isotropy arise at the level of 1 part in 10^; these 
fluctuations arise because of the non-uniform distribution of matter at (and after) the 
surface of last scattering. 

It is therefore not unreasonable to assume that the stochastic gravity wave background 
is also isotropic. However this assumption may not be true. For example, suppose that 
the dominant source of stochastic gravity wave background is a large number of unre- 
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Figure 2. A conformal diagram showing a spatially-flat Friedman- Robertson- Walker cosmo- 
logical model. The past light cone of a present-day observer (Here, Now) intersects a surface at 
t = 10~^^ seconds in a large 2-sphere, of which two points (a,b) are shown. The future light cone 
of the initial singularity intersects the same surface in a much smaller 2-sphere, of which two 
points (a,c) are also shown. The number of independent, uncorrelated horizon volumes A'horizon 
which contribute to the gravitational radiation arriving at a detector today is given by the ratio 
of the areas of the larger 2-sphere to the smaller one. 



solved white dwarf binaries within our own galaxy. Because our galaxy is bar or spiral 
shaped (and not spherical) if we assume that the white dwarf binaries are distributed 
in space in the same way as the matter in the galaxy, then the stochastic background 
will have a distinctly anisotropic distribution, and will form a "band in the sky" dis- 
tributed roughly in the same way as the milky way. It is also possible for a stochastic 
gravity wave background of cosmological origin to be quite anisotropic. Recent work 



(Allen fc Ottewill 1996) has shown how data from the LIGO detectors could be ana- 



lyzed to search for such anisotropics. Of course in this case, one is left with the problem 
of explaining why the CMBR is isotropic, but the gravity wave background is not. It is 
possible to conceive of such mechanisms, but within the scope of these lectures, we will 
not consider their effects. 

The assumption that the stochastic background is stationary is almost certainly jus- 
tified. Technically this means that the rt-point correlation functions of the gravitational 
wave fields depend only upon the differences between the times, and not on the choice 
of the time origin. Because the age of the universe is 20 orders of magnitude larger than 
the characteristic period of the waves that LIGO, VIRGO, and the other facilities can 
detect, and 9 orders of magnitude larger than the longest realistic observation times, it 
seems very unlikely that the stochastic background has statistical properties that vary 
over either of these time-scales. 

The final assumption we make is that the fields are Gaussian, by which we mean that 
the joint density function is a multivariate normal distribution. For many early-universe 
processes which give rise to a stochastic background, this is a reasonable assumption, 
which can be justified with the central limit theorem. Let's sketch the argument out. 

We will show in Section 2.4 that if the stochastic background arises from processes that 
take place in the early universe, the characteristic time (proper time after the big bang) 
at which the gravitational radiation was emitted was about t = 10~^^ sec. Consider 
the conformal diagram (Hawking & Ellis 1973), shown in Fig. |^. As Fig. || graphically 
illustrates, a detector today located at the spacetime point labeled "Here, Now" observes 
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radiation produced aXt — 10^^^ seconds after the big bang by an extremely large number 
-^'horizon of independent horizon volumes. We can estimate the number of these horizon- 
sized volumes on the surface t = 10^^^ seconds as follows. 

Let's assume that the universe is fc = (spatially flat) and radiation dominated from 
the time that the gravitons were produced (ti — 10~^^ sec) until the present time (to = 
10^^ sec). This is a reasonable assumption for this type of calculation; although the 
universe did recently become matter dominated (at a redshift Z of a few thousand) this 
has only a small effect on the final answer. The redshift Z of the t = ti constant-time 
surface is then given by 1 + Z = (10^^/10^^^)^/^ ~ 10^". Let's work in a "conformal 
time" coordinates where ds^ = a^(r])(—dr]^ + dx^). For a radiation-dominated universe 
a(?7) — rj, and we can take the present time to have rj ^ rjo and the time at which 
gravitons were produced to have 77 = 771 . We begin by considering the intersection of the 
forward light cone of the initial singularity with the surface rj — rji. This intersection 
forms the small two-sphere denoted by the points "a" and "c" in Fig. |[ this two-sphere 
has area Agmaii = 47r?7ja^(?7i). In similar fashion, the past light cone of the present-day 
detector intersects the surface 77 = 771 in the large two sphere denoted by the points "a" 
and "b", which has area Abig = 47r(r/o — 771)^0^(771). The ratio of these two areas is the 
number of horizon volumes visible to the detector over the entire sky; it is 



^big _ (770 - 7 71^^ 



iV,,,,_ = ^ ^ ^Z'^ 10^^ (2.10) 



If we assume that the processes producing gravitational waves from each of these separate 
horizon volumes act independently, then it follows immediately from the central limit 
theorem that the amplitude of the radiation arriving at the detector, which is the sum of 
the amplitudes of the radiation produced by each of the separate volumes, is Gaussian. 

2.4. When is a stochastic background produced? 

Suppose that the LIGO or VIRGO detectors do detect a stochastic background of grav- 
itational radiation. From what epoch, in the history of the universe, does this radiation 
date? 

To give a general answer to this question, we need to make some assumptions about the 
universe. The most reasonable thing to do is to adopt the "standard model" of cosmology. 
This cosmological model consists of a spatially-flat (fc = Friedman-Robertson- Walker) 
cosmological model, in which the equation of state is dominated by massless (or highly 
relativistic) matter for redshifts greater than Z^q « 6000, and is dominated by massive 
pressureless particles (dust!) for smaller redshifts. 

The question posed above has two answers, one mundane, the other quite exciting. 
The mundane possibility is that the radiation was produced at or near the present epoch 
(say, within a redshift Z < 4) hy many unresolved separate sources, such as white dwarf 
binaries or supernovae. For example, in the case of the electromagnetic background 
radiation, one finds that below about 300 MHz the spectrum is dominated by (recent) 



emission from our own galaxy (Weinberg 1972). The more exciting possibility is that 



such gravitational radiation is a cosmological fossil. In the case of the electromagnetic 
spectrum, at frequencies above 1 GHz, this is indeed the case. 

Note that we must not discount the possibility that the stochastic background in the 
LIGO frequency band was produced near the present epoch. Recent work (Blair fc Lu 1996)] 



has shown that if the efficiency with which supernovae convert their mass into gravita- 
tional radiation exceeds 10~^ then the cosmological distribution of such supernovae could 
be a dominant (and detectable) source of stochastic background in the LIGO frequency 
band. 
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Figure 3. This graph shows the two lengths, as functions of cosmological time (both axes 
are logjg). Today, we are at the far right of the graph, t ^ 10^'' seconds after the big bang. 
The solid curve shows the physical wavelength of a wave that LIGO or VIRGO might detect, at 
100 Hz. Today this wavelength is about 3 x 10* cm; in the past this wavelength shrinks because 
of redshifting, in proportion to the cosmological scale factor, first as t^^^ when the universe is 
matter dominated, then as t^^^ when the universe is radiation dominated. The dotted curve 
shows the characteristic spatial size (and age) of the universe, the Hubble length, which today 
is r; 10^* cm. This function grows cx t. The intersection of these two curves determines the time 
at which the primordial production of gravitons in the LIGO/ VIRGO band took place. 



Let us assume that the LIGO or VIRGO detectors do detect a stochastic background 
of gravitational radiation, and let us also assume that this radiation has a cosmological 
(a long time ago) rather than a local (recently) origin. The time at which this radiation 
was produced can be determined by simple physical argument. The production of grav- 
itational radiation is no different than the production of any other type of radiation. If 
we want to produce electromagnetic radiation at 1 KHz, we need to take electric charges 
and vibrate them at 1 KHz. The same holds for gravitational radiation; waves of a 
certain frequency are produced when the characteristic time for the matter and energy 
in the universe to shift about is comparable to the period of the waves. The graph in 
Fig. P shows two length scales, as functions of cosmic time. As can be seen from the 
diagram, at about 10~^^ seconds after the big bang the characteristic expansion time 
(the Hubble time) was about comparable to the period of a wave, which when redshifted 
to the present epoch, falls within the bandwidth of LIGO /VIRGO sensitivity. 

The most exciting prospect in the study of the stochastic background is the possibility 
of detecting radiation produced at this very early epoch of t « 10~^^ sec. To give a sense 
for the scales, a 100 Hz graviton detected by the LIGO experiment was produced at that 
early epoch with an energy of order 10 MeV. The mean temperature in the universe at 
that time was of order lO"^ GeV. For this reason, detecting a background of cosmological 
origin would give us a glimpse of the universe at much earlier times than we can obtain 
in other ways. 

Comparable calculations can be done for the proposed Laser Interferometer Space 



Antenna (LISA) experiment (Hough et al. 1995), which is designed to detect waves with 
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frequencies in the range « 10^^ Hz to « 10~^ Hz. In this case, one finds that waves of 
cosmological origin with present-day frequency « 10~^Hz were produced at time of order 
10"^"* sec after the big bang; at the time of their production the energy of the gravitons 
was « 0.1 eV. The overall temperature of the universe at this time was k, 10^ GeV - a 
particularly interesting epoch, because the electro- weak phase transition ocurred around 
that energy scale. 



3. Experimental detection methods and their sensitivity 

3.1. How to detect a stochastic background (crude version) 

In these lectures, I'd like to answer the question "how does one actually use the output 
from the gravitational wave detectors to find (or place upper limits on) a stochastic 
gravity wave background." Later, in Section 3.3 we will re-visit this question and give a 
very precise (in fact, optimal) answer to this question. However my goal here is to give 
a rough understanding of the method and its limitations. 

Now in general, there are two possibilities. First, when a detector is operating, its 
intrinsic noise level might be lower than the random stochastic signal arising from the 
gravitational wave background. In this case, no clever signal detection strategies are 
required. One simply "turns on" the detector, studies its output for a few seconds, and 
concludes that a stochastic background is present. In a crude sense, this is how Penzias 
and Wilson discovered the electromagnetic CMBR; the signal it produced was greater 
than the intrinsic noise of their detector. Of course it was not quite so easy - Penzias 
and Wilson had to spend an entire summer studying the intrinsic noise in their radio 
telescope before they could be sure that the excess noise they were observing was of 
cosmological origin! 

Let us assume that the signal levels produced by the stochastic background of grav- 
itational radiation are not this large, so that the output of a detector is dominated by 
the noise produced within the detector rather than by the signal due to the stochastic 
background itself. In this case, as was first shown by Michelson (1987) one can still de- 
tect the signal due to a stochastic background, by correlating the outputs of two different 
detectors. The principal requirement that must be met by the two detectors is that they 
must have independent noise. Later in these lectures we will be more specific about what 
is required. 

In the case where we have two detectors, one can isolate the signal due to the stochastic 
background by correlating their outputs. As can be seen from Fig. ^, a gravitational plane 
wave passing by the earth excites a pair of detectors in coincidence when the positive 
or negative amplitude part of the wave is passing by both detectors at the same time 
(in the rest frame of the detectors). It excites the detectors in anti-coincidence during 
moments when when the positive amplitude part of the wave is passing by one detector, 
and the negative amplitude part of the wave is passing by the other detector, or vice 
versa. Provided that the wavelength of the wave is much larger than twice the distance 
between the two detectors, on the average the detectors will be driven in coincidence. 
In the case of the LIGO detector pair, below a frequency of 64 Hz, the Hanford and 
Livingston detectors will be driven (on the average) in coincidence by an unpolarized 
isotropic stochastic background. This can be seen by inspecting Fig. ||; the first zero 
of this function lies at / = 64 Hz. We note that this result is slightly higher than the 
frequency (50 Hz) associated with a wavelength of twice the detector separation (6000 
Km) because the overlap reduction function is a sum of Bessel functions (B.ll) which 
does not vanish at 50 Hz. 
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Figure 4. Shown is a the surface of the earth 15° < latitude < 75°, and 

— 130° < longitude < 20°, including the LIGO detectors in Hanford, WA (LI) and Livingston, 
LA (L2) and the locations of the VIRGO detector (V) in Pisa, Italy and the GEO-600 (G) 
detector in Hannover, Germany. The perpendicular arms of the LIGO detectors are also illus- 
trated (though not to scale). A plane gravitational wave passing by the earth is indicated by 
successive minima (troughs), zeros, and maxima (peaks) of the wave. As this wave passes by the 
pair of LIGO detectors, it excites the detectors in coincidence at the moment shown, because 
both detectors are driven negative by the wave. During the time when the zero is between LI 
and L2, the two detectors respond in anti-coincidence. Provided that the wavelength A of the 
wave is larger than about twice the separation distance (2998 Km) between the detectors, on 
the average they are driven in coincidence more of the time than in anti-coincidence. 



The basic idea (the optimal implementation is elaborated in Section 3.3) is to multiply 
together the outputs of the two independent detectors, and integrate the result. For 
the purpose of this section, imagine that two detectors are located at the same point, 
and have identical orientations. In our case, for frequencies well below 64 Hz, the LIGO 
detectors may be thought of as "being at the same point" . In similar fashion, although 
the LIGO detectors are not identically oriented (they are 27.2° apart as seen from the 
center of the earth, and thus can not lie in the same plane) their arms are as close to 
parallel as is possible given their separation. 

In this simplified case the output of the first detector is 

si{t)^hi{t)+n^{t), (3.1) 

where hi is the strain due to the stochastic background and ni is the intrinsic noise of 
the first detector. In similar fashion, the output of the second detector is 

S2{t)^h2{t)+n2{t). (3.2) 

Since we are assuming that the two detectors have identical locations and orientations 
then the gravity- wave strains are identical: hiit) — h2{t). We can form a "signal" by 
multiplying together the outputs of the two detectors and integrating: 

.T/2 

S^<si,S2>= si{t)s2{t)dt. (3.3) 



-T/2 

Here T is the integration time (realistically, a few months). There are two different 
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cases to consider. As we explained earlier, the case where the detector noise Ui is small 
compared to the signal Si not very interesting, as one can find the background easily in 
that case. So we'll assume that the signal is small compared to the detector noise, in 
which case we can write 

S = < hi,h2 > + < ni,h2 > + < hi,n2 > + < ni,n2 > (3.4) 
« < /ii, /i2 > + < f^2 > • (3.5) 

where we have dropped terms like < ni,/i2 > that are smaller than < ni,n2 > but 
statistically identical, because ni and /12 are uncorrelated in the same way as ni and 
712. (To simplify, we assume that there are no correlated sources of noise such as seismic 
activity) . 

The two terms appearing in the signal S grow in different ways as the observation time 
T increases. The first term < hi,h2 > increases linearly with time T, because hi — /12. 
The second term < ni,n2 > can be thought of as a random walk on a one-dimensional 



line, and so on the average it grows oc \/T. We will show in Section 3.3 that 



</ii,/i2 >« |/i(/)|2A/T(xr!(/)A/r (3.6) 
<ni,n2 >« |n(/)|2yA7r (3.7) 

where / is the central frequency at which the detector is sensitive. A/ is the effective 
bandwidth, and Q is the Fourier transform of Q. Because the "signal" < /ii, /i2 > grows 
like r, and the noise < rti,rt2 > only grows as Ti/^ with enough observation time one 
can in principle detect a gravitational wave stochastic background buried in any level of 
detector noise. Setting the "signal" term greater than the noise term gives a minimum 
detectable level of VI: 

^ ^minimum detectable ^ / A -f ^ ' \ / 

Thus, even if the detector noise is much larger than the expected stochastic gravity wave 
strain, by making a long enough observation, one can in principle observe even very low 
levels of stochastic background. 

Let me stress one more time that this crude analysis is intended pri mar ily to present 
the concepts and basic ideas, and will be treated rigorously in Section |3.3|. 



3.2. The overlap reduction function j{f) 

To provide a rigorous treatment of the signal analysis, we must take into account the 
reduction in sensitivity that results from two effects (1) the non-parallel alignment of the 
arms, and (2) the separation time delay between the two detectors. These two effects 
reduce the sensitivity of a stochastic background search; they mean that hi and /12 are 
no longer equal; the overlap between the gravity wave strains in the two detectors is only 
partial. 

To quantify these effects we introduce the overlap reduction function "f{f) first calcu- 
lated in closed form by Flanagan (1993). This is a dimensionless function of frequency /, 
which depends entirely on the relative positions and orientations of a pair of detectors. 
The overlap reduction function 7(/) equals unity for co-aligned and coincident detectors, 
and decreases below unity when the detectors are rotated out of co-alignment (so they 
are sensitive to different polarizations) or shifted apart (so that there is a phase-shift 
between the signals in the two detectors). Later in these lectures (in Section |3.3| ) we will 
see how this function comes up naturally in the determination of signal-to-noise ratios. 
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Figure 5. The overlap reduction function 7(/) for the two LIGO detector sites. (The 
horizontal axis of the left-hand graph is linear, while that of the right-hand graph is logio.) 
The overlap reduction function shows how the correlation of the detector pair to an unpolarized 
stochastic background falls off with frequency. The overlap reduction function has its first zero 
at 64 Hz, as explained earlier. It falls off rapidly at higher frequencies. 



The overlap reduction function is defined by [(Flanagan 1993) 

5 



lif) 



(3.9) 



Here is a unit-length vector on the two-sphere, Aa' is the separation between the two 
detector sites, and F^'^ is the response of detector i to the -f or x polarization. For the 
first detector {i — 1) one has 

1 



(3.10) 



where the directions of a single detector's arms are defined by Xi and Yf-, and e^,l^ (Cl) 
are the spin-two polarization tensors for the "plus" and "cross" polarizations respectively. 
These quantities are defined in equations (A_7,A_8) of Appendix 



It is not hard to see how the overlap reduction function j{f) arises. As we show in 
detail in Appendix ^ 7(/) is proportional to the averaged product of the strains at two 
different detectors, when those detectors are driven by an isotropic unpolarized stochastic 
background of gravitational radiation at frequency /. This radiation is arriving at the 
detectors uniformly from all points on the celestial sphere, hence the integral over the 
direction {l of the wave- vector of the arriving radiation. 

Let us look at each term in ( |3.9[ ). (1) The overall normalization factor is chosen so 
that if I Ax I — (or equivalently, / = 0) and the two detectors are coincident, and the 
two detectors have parallel pairs of arms, then 7^1. (2) The exponential phase factor 
is the phase shift arising from the time delay between the two detectors for radiation 
arriving along the direction tl. (3) The quantity F^ is the response of the first detector 
to the "-I-" polarization wave, and similarly for detector 2 and the "cross" polarization. 



In Section 3.3 this overlap reduction function will appear after a rigorous calculation; 
these remarks are meant primarily for motivation. 

A graph of "f{f) is shown in Fig. ^ for the two LIGO sites in Hanford, WA and 
Livingston LA. The overlap reduction function is negative as / — s- 0. This is because, 
comparing the two LIGO sites, the arm orientations are not parallel but are rotated 
relative to one another by 90°. This can be seen immediately in Fig. ^. If the Livingston 
detector arms (denoted L2) were rotated by 90° in the clockwise direction, this would 
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change the overall sign of 7(/), but otherwise leave its dependence on / unchanged. The 
magnitude of 7(0) is not unity, because the plane of the Hanford detector and that of 
the Livingston detector are different. Thus the arms of the two detectors are not exactly 
parallel, and hence |7(0)| = 0.89 rather than 1. 

The precise form of the overlap reduction function is not important here. However it 
can be calculated in closed form (Flanagan 1993), and may be expressed as the sum of 



three Bessel functions: 

l{f) = Ajo[a) + B hC — 2- 3.11 

where a — 27r/|Aa;|/c depends upon the detector separation and the frequency, and 
A, B, C are constants whi ch only depend u pon the locations of the detectors and the 
orientations of their arms [(Flanagan 1993) . For example, for the two LIGO detectors, 



one finds A = -0.124842, B = -2.90014, and C = 3.00837. 

3.3. Optimal filtering 



I now return to the subject that we addressed in Section 3.1: how do we detect a stochastic 
background? In this section, we work out the precise answer to this question. This type of 
analysis was first carried out by Michelson (1987) then improved by Christensen (1992) 
and further refined by Flanagan (1993). The analysis given here appears in Bayesian 
form in the work of Flanagan (1993); the treatment given here is the somewhat simpler 
"frequentist" version. The two different approaches yield exactly the same final results, 
because the signals are Gaussian. 

Let us write the most general possible form of the correlation "signal" between two 
detectors in the form 

.T/2 .T/2 

S= dt dt'si{t)s2{t')Q[t-t'). (3.12) 

J-T/2 J-T/2 

The strains Si and S2 are the signal streams from (say) the Hanford and Livingston 



LIGO detectors, defined by (3.1-3.2), and Q{t ~ t') is a real filter function. Because we 
are assuming that the statistical properties of our signals and noise are stationary, the 
"best choice" of filter function Q can only depend upon the difference At — t — t' oi the 
two times. Our goal, in this section, is to find the optimal choice of filter function Q(At) 
in a rigorous way. 

The optimal choice of filter function Q{At) depends upon the locations and orienta- 
tions of the detectors, as well as on the spectrum of the stochastic background, and the 
noise characteristics of the detectors. If the two sites are very close together (compared 
to the characteristic wavelength at which the detectors are sensitive) , have identical ori- 
entations, and the detectors have a narrow bandwidth, then the optimal choice of filter is 
a Dirac delta function Q{t — t') — 5(t — t'). In this case, we just multiply the signals and 



integrate, precisely as discussed earlier in Section Even if none of these assumptions 
holds, you will see that the optimal choice of filter function Q falls rapidly to zero for 
time delays At whose magnitude is large compared to the light travel time d/c = 10^^ 
seconds between the two sites. We are thus justified in changing the limits above, to 
obtain 

/T/2 .00 
dt / dt'si{t)s2{t')Q{t-t'). (3.13) 
-T/2 J-co 

Since the observation time T will be at least several months, one can see immediately 
that the signal S gets most of its support in the integral ( 3.12| ) from the region near 



the diagonal line t = t' in the region [-r/2,r/2] x [-r/2,T/2]. Hence the value of S 
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grows linearly with integration time T, in perfect agreement with the earlier "intuitive" 
explanation in Section 



3.1 



The reader should note that with our definition, the "signal" 5 is a quantity quadratic 
in the strain h\ this is not the same as the convention frequently followed in analyzing 
signals from chirp or periodic sources, where the "signal" is generally taken to be a 
quantity linear in the strain h. For this reason, the expressions that we will derive for 
quantities such as the signal to noise ratio are in some ways analogous to the squares of 
the corresponding quantities for chirp or periodic signals. 



We can re- write (3.13) in the frequency domain. Our convention for Fourier transforms 

is 

/•oo 

.g(/)=/ dte-^^^f'git). (3.14) 



Here, and elsewhere in these lectures, * denotes complex conjugation, and within this 
section one may assume that Si{t) vanishes at some very early time (say, before the 
detector was built) and at some very late time (say, after the detector is decommissioned) 
so that the Fourier transforms are well-defined. Because we have assumed that the filter 
function is real, Q{—f) = Q*{f)- Identical reality conditions hold for the signals Si, the 
strains /i^, and the detector noise fii. Using the Fourier transform, we may write the 



signal ( |3.13| ) as 

5=/ df I df StU - nsUf)S2{f')Q{f). (3.15) 



The function St which appears on the r.h.s. of (3.15) is a "finite time" approximation 
to the Dirac delta function, defined by 

Srif)^ r^te'^-f^^'-^^^, (3.16) 

J-T/2 '^J 

which reduces to the Dirac 5-function 5{f) in the limit T — > oo. Note however that for 
finite observation times T one has 6t{0) — T. We now consider the stochastic properties 
of the signal 5*. In order to find the optimal choice of the filter function, some additional 
assumptions are needed. 

As discussed in detail earlier, we will assume that the stochastic background is isotropic, 
unpolarized, and Gaussian. Under these assumptions, the expectation value of the 
Fourier amplitudes of the strain signal and the strain noise are derived in Appendix ^ 

ChUmif)) = 6if f)^\f\-^n,^{\f\)^{\f\), and (3.17) 

{nKmAf')) = ls{f /')%P.(I/I). (3.18) 

Here the indices i,j label the sites, so for example i = 1 is Hanford Washington and 
z = 2 is Livingston Louisiana. The power spectral noise density function Pi{f) which is 



defined by (3.18) is a real non-negative function, defined with a factor of 1/2 to agree 
with the standard definition used by instrument-builders, so that the total noise power 
is the integral of P{f) over all frequencies from to oo (rather than starting at — oo). 
You will notice that we have assumed that the detector noise at the two different sites 
is completely uncorrelated. This assumption significantly simplifies the analysis, but is 
not entirely realistic (Christensen 1992; Flanagan 1993). 

The predicted detector noise for the initial and advanced LIGO is shown in Fig. O. 



This figure is taken from the LIGO design goals (Abramovici et al. 1992), Because the 



4 kilometer detectors are identical, they have the same power spectral noise densities. 
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Log Freq, Hz 

Figure 6. The predicted noise power spectra of the initial and advanced LIGO detectors. 
The horizontal axis is logio of frequency /, in Hz. The vertical axis shows logio(P(f)/sec)^/'^, 
or strain per root Hz. These noise power spectra are the published design goals. The bumps 
appearing in the low-frequency part of the advanced LIGO noise curve are obtained by folding 
measured seismic noise data with the predicted transfer function of the seismic isolation (stack) 
system. 



P{f) — Pi{f) — P2{f)- There are at least three different sources of detector noise which 
dominate at different frequencies. At the low frequency end, the dominant source of 
noise is seismic noise. In the advanced version of LIGO, this noise will be reduced by 
improving the mechanical filters that isolate the detector from the ground. In the mid- 
band region, the noise is primarily due to thermal motion of the masses and mirrors, 
which will be reduced by using improved materials and design. In the high-frequency 
part of the spectrum, the noise is dominated by photon shot noise, which can be reduced 
by increasing the laser power and improving the optical components. 

Although in the case of the LIGO detectors, the noise power spectra of the detectors 
are expected to be the same, in general, we are interested in correlating the outputs of 
detectors with different designs (say the LIGO and VIRGO detectors). So we will keep 
our treatment fairly general, and treat Pi{f) and P2{f) as distinct functions. 

In order to find the optimal choice of filter function Q, we need to choose the quantity 
to maximize. The natural choice here is to maximize the ratio of signal-to-noise. The 



expected value of the signal is obtained by taking the expectation value of (3.15), making 



use of (3.17). Because the noise in each detector is uncorrelated with the other detector, 



and is uncorrelated with the gravity-wave strain ft,, one obtains 

/•oo />oo 

{s)^ df df' 5t{i ~ nKifYMnmn (3.19) 



—00 J —00 

2 /"OO 



^^'Ti d/7(l/l)l/r'f^gw(|/|)Q(/). (3.20) 



207r 



2 ^ 



—00 



The factor of T on the r.h.s. arises from evaluating (5t(0). To determine the amount of 
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noise N in the detector, we consider the variation of the signal S away from its mean 
value: 



N = S-{S). 



(3.21) 



The quantities on the r.h.s. are easily determined, if we assume as in (3.4) that the noise 
ni{t) each detector is much larger in magnitude than the gravitational strain hi(t). In 
this case, one may simply r eplac e the detector output signal (gravity strain + noise) by 
the detector noise alone in ( 3.15 ), obtaining 

/■oc 



N ' 



df 



df StU - f)nl{f)Mf')Q{n 



(3.22) 



This approximation is consistent with the definition of N; a brief calculation shows that 
the expected value of N vanishes, because the noise in the two detectors is assumed to 
be uncorrelated, 

(N) = 0, (3.23) 

a result which also follows directly from the definition of N. However the r.m.s. value of 
the expected noise is not zero, because if we take the expectation value of the square of 
(3.22) and use ( 3.18D one obtains 

(N^) = (^2^ - (5)2 

1 /"^ 



df 5Uf - npiiifmnri) wo 



(3.24) 
(3.25) 



For long observation times (say T > 1 sec) the finite-time delta function (5t(/ — /') is 
sharply peaked over a region in / — /' whose size « 1/T is very small compared to 
the scale on which the functions Pi, P2,Q are changing. In this case we are justified in 
replacing one of the finite-time delta functions above by a Dirac delta function, obtaining 

rj-i />00 

{N') = - / rf/Pi(|/|)P2(|./|) Q{.f) 



(3.26) 

Notice that the expected r.m.s. noise {N'^)i grows as VT, in complete agreement with 
the simple intuitive explanation given earlier in Section 3.1 



The problem now is to find the choice of filter function Q which maximizes the signal 
to noise ratio 



(7V2 



(3.27) 



This turns out to be remarkably straightforward, if we introduce an inner product be- 
tween two functions of frequency, A{f) and B{f). This inner product takes as inputs 
two arbitrary complex functions A{f) and B{f); the output {A, B) is a single complex 
number. We define this inner product by 



(AB) 



dfA*{f)Bif)Pii\f\)P2i\f\)- 



(3.28) 



It is straightforward to see that this inner product is positive definite; since < Pi{f) 
the norm of a function [A, A) vanishes if and only if A{f) vanishes. In terms of this inner 
product, we can express the expected signal and noise as 



_ 7(l/l)»gw(|/|) 3Hl 
^^"^^'l/PA(|/|)P2(|/|)^ 207r2'' 



(3.29) 
(3.30) 




Figure 7. Optimal filter functions Q{f) for the initial and advanced LIGO detectors, for 
f2gw = constant. In this case Q{f) is proportional to jTp^jj^- The initial LIGO filter is 
normalized to have maximum magnitude of unity; the advanced LIGO filter is shown with an 
identical scale. It is much larger in magnitude because the power spectral noise density P{f) of 
the advanced detectors is much smaller than in the initial detectors. 



So the problem is to choose the function Q which gives the largest possible value of 

C\2 /c\2 /Qrr2\2 ( f) 7(l/l)»gw(l/l) \2 

n) -{N^)-[lO.^) iQ,Q) ■ ^'-''^ 

But this is trivial! Suppose you were given a fixed three-vector A, and asked to find the 
three- vector Q which maximized the ratio '■^'^^ . Since this is just the squared cosine 

of the angle between the two vectors, the ratio is obviously maximized by choosing A 
and Q to point in the same direction! The problem of maximizing (B.31) is identical; the 
solution is to choose 

^/ _ 7(|/I)^gw(|/|) „„x 

\f\'Pii\f\)P2i\f\y ^ ^ ^ 

Of course the overall constant appearing in the normalization of Q has no significance, 
since it does not affect the signal-to- noise ratio ( ^.31 ). We have now found the form of 
the optimal filter for a stochastic background search. 

These optimal filters are shown in Fig. ^ for both the initial and advanced LIGO 
detectors, for a flat spectrum rig„ — constant. Effectively, the filter includes into the 
signal only those frequencies at which the detector is most sensitive, weighting each 
frequency by a positive or negative factor which reflects the relative phasing (at that 
frequency) of the signals at the two detectors. 

Having now found the optimal choice of filter function, it is straightforward to calculate 
the signal-to- noise ratios for a given pair of detectors. Substituting the optimal choice of 
Q given by (3.32) into (3.31) gives the optimal signal-to-noise ratio: 

'"'Tr<if:y]y!i (^a) 



NJ 50»« Jo f'PiU)PM) 

Thus the signal-to-noise ratio of a correlation experiment to detect a stochastic gravita- 
tional wave background is determined by the following quantities. 

iJo — Hubble expansion rate, sec~^. 

T = Integration time, sec. 

li.f) = Overlap reduction function, dimensionless. 
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rigw = Fractional energy density in gravitational waves. 
Pi = Noise power spectral density of detector i (one sided), sec. 

One of the curious things about this formula for the signal-to-noise ratio is that it depends 
upon the spectrum of gravitational waves fig^; a function which we do not know! What 
this means in practice is that rather than having a single "optimal filter" one needs to have 
a set of such filters. The bandwidth of the experiment is in a good sense defined by the 



factor (f )P2 (/) appearing in the integrand of (3.33); the bandwidth of the experiment 
is the range of frequencies over which this factor is large (and hence the noise is small). 
The signal-to-noise ratio gets its significant contributions from within this bandwidth. 
Within this range of frequencies, which will not be very wide for the first generation of 
detectors, it is reasonable to assume that the form of the figw is a power-law, proportional 
to Thus, one can construct a reasonable set of optimal filters by choosing a number of 
different values of a such as a = —4, —7/2, • • • , 7/2, 4 and analyzing the signal separately 
for each filter. The filter for a = is shown in Fig. 0. 

In the next section, we will look at the values of this signal-to-noise ratio for the 
forthcoming detectors, and estimate the detectable levels of figw. 

3.4. Signal to noise ratios, and the minimum detectable figw 
One can see from ( 3.15| ) that in a measurement over the time interval T the signal S is 



a sum (over / and /') of many independent random variables (the Fourier amplitudes of 
the signals). This is because the amplitudes si(/) and si(/') are only correlated when 



1/ — /'I < 1/T « 10~^ Hz, and the bandwidth over which the integral in ( 3.15 ) gets its 
major contribution is ~ 100 Hz wide. Thus S is the sum of ~ 10^ random variables, and 
the value of S" in a set of measurements over independent time intervals is a Gaussian 
random variable. If the signal S is measured many times, each time over a distinct period 
of length T, the values of S will have a Gaussian normal distribution. The mean value of 
this distribution is (5*); the width or standard deviation of this distribution away from 
its mean value is (TV^)^ . 

In order to reliably detect a stochastic background, we need to be able to state, with 
some certainty, that the value of Ogw is greater than zero. Equivalently, one needs to 
be able to state with some desired reliability that the observed positive mean value of 
S could not have resulted from detector noise but instead must have resulted from a 
stochastic background. 

One way of doing this is to compute the probability that a random variable with a 
Gaussian normal distribution will lie within a given range. For example, the probability 
that a Gaussian random variable 5* with mean {S) and variance = {N"^) will he within 
the range [Si , Su] is 

1 /"^" f — ix— < S >)'^\ 

Probability = , / dx exp ^ — — . (3.34) 

This integral can be expressed in closed form in terms of the error function Erf. In this 
way, one can easily show that 



with 68% probability S £ [{S) - y/W), (S) + y/{N^)] 
with 90% probability S G [{S) - 1.65^^^, (5) + 1.65^(^^)1 
with 95% probability S £ [{S) - 2y/JJP), (S) + 2^/(JP)] 
with 99.7% probability S G [{S) - 3\/Jn^, {S) + 3y/JTP)]. 
Thus, to detect a stochastic background with > 90% confidence, we need to ensure that 
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Figure 8. The expected "signal" as a function of frequency for a flat spectrum Jlgw — constant 
in the initial and advanced LIGO detectors. The function displayed here is the integrand of Eqn. 

3.20 ; it is proportional to -j^-pf^j^- The initial LIGO signal is normalized to have maximum 
magnitude of unity; the advanced LIGO signal is shown with an identical scale. It is much 



the advanced LIGO signal is shown with an identical scale, 
larger in magnitude because the power spectral noise density of the advanced detectors is much 
smaller than in the initial detectors. 



the signal-to- noise ratio (-^) > 1.65, so that there is less than a ten percent probability 
that S is less than 0. 

As an example, let's analyze the sensitivity of the first and second generation of LIGO 
detectors, to a flat spectrum ilg^ c>c /°. The design goals (Abramovici et al. 1992) for 
the noise power spectral densities of the detectors are shown in Fig. |^. The minimum 
flat spectral density detectable in T = 4 months of integration, with 90% confidence, is 
given by 



90% 
gw 



[ifide 



9m 



507r4 



2.8 X IQ-^h^"^ 



fPiif)P2if) 
iQQ for initial LIGO 



1.65 ^Pi-T / df ^f,J,y,l, (3.35) 

(3.36) 

= 2.8 X lO^^/ij^oo for advanced LIGO. (3.37) 

We will see in the coming lectures that this is a very useful level of sensitivity, particularly 
for the advanced detectors. 

One might wonder why the upgrade from the initial LIGO detectors to the advanced 
detectors results in such a large increase in sensitivity. The reason can be easily seen 
by comparing the optimal filter functions Q{f) for the initial and adv anced detectors 
shown in Fig. 0, and the expected contributions to (5) given by ( 3.20 ), as a function 
of frequency, shown in Fig. |[ It is easy to see that because the advanced version of 
the LIGO detector will have much less low-frequency noise than the initial detector, the 
expected signal in the LIGO detector lies within the < / < 64 Hz region where the 
overlap reduction function "f{f) has the greatest magnitude. 



4. What is known (experimental/observational facts) about the 
stochastic background? 

So far in these lectures, we have discussed notation which may be used to describe 
the spectrum of a stochastic gravitational wave background, and some of the expected 
properties that such a background might have. We also discussed, in detail, the optimal 
way to correlate the output of two interferometers to search for such a background. 
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In this section, I'd like to show you what is actually known about the stochastic 
background: the observational facts. As you will see, there are strict limits on the 
stochastic background in only a couple of frequency ranges, but other than that, only 
one very general constraint. In particular, there are really no strong constraints on O 
in the frequency range of interest for ground based interferometric experiments or bar 
detectors (1 Hz < / < 10"^ Hz) or space-based experiments (10^"* Hz < / < 10^^ Hz). 

4.1. CMBR isotropy and the gravitational stochastic background 

The strongest constraint on Q comes from the high degree of isotropy observed in the 
Cosmic Microwave Background Radiation (CMBR). In particular, t he one-year (Smoot et 



al. 1992; Bennett et al. 1992; Wright 1992) and two-year data-sets [(Bennett et al. 1994) 
from the Cosmic Background Explorer (COBE) satellite place very strong restrictions 
on f2 at very low frequencies. The four-year data has also appeared, although it is not 



yet in print [(Bennett et al. 1996) 



The temperature fluctuation ST of the CMBR away from its mean value of T = 2.73 K 
varies from point to point on the sky, and may be expanded in harmonic functions on 
the celestial sphere: 

^(ii) = ^a,,™y,™(f]). (4.1) 

Im 

Here f2 denotes a point on the two-sphere, and the ai^m are the multipole moments of 
the CMBR. Since by definition the mean value of the l.h.s. vanishes, the sum begins at 
/ = 1: 

OO / 

E-EE- (4-2) 

Im l — l m— — l 

The COBE data permitted a determination of the rms quadrupole moment of these 
fluctuations, 



= T 



\ 



^ %C = 14.3«JmK. (4.3) 

These in turn place rather tight constraints on the spectrum of gravitational radiation 
at very long wavelengths. 

The connection between the temperature fluctuations of the CMBR and gravitational 
waves arises through the Sachs Wolfe effect. Sachs and Wolfe (1967) showed how varia- 
tions in the density of the cosmological fluid and gravitational wave perturbations result 
in CMBR temperature fluctuations, even if the surface of last scattering was perfectly 
uniform in temperature. In particular, they showed that to lowest order the temperature 
perturbation at a particular point on the sky is given by 

^(f)) = U dX^Kr. (4.4) 

^ J null geodesic 

Here the integral is taken over a null geodesic path which leaves our current spacetime 
point heading off in the spatial direction defined by J7 and going back to the surface 
of last scattering of the CMBR, and A is a particular choice of affine parameter along 
that null geodesic. In conformal coordinates the metric perturbation is given by Sgij = 
a'^{rj)hij, and r is an outwards radial spatial coordinate from our current spacetime 
point. The effect of a long wavelength gravitational wave is to shift the temperature 
distribution of the CMBR on the celestial sphere away from perfect isotropy. The fact 
that the temperature perturbations are quite small (a few parts in 10^) implies that 
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Table 1. Two of the millisecond pulsars studied by Taylor's group, and then used to place 

limits on the stochastic gravitational wave background. 

Pulsar Name Time Observed t Period Timing Residuals At 
PSR B1855+09 7 Years 5.4 msec no residuals 

PSR B1937+21 8 Years 1.6 msec few ^sec (total) 



the perturbations caused by the gravitational waves can not be too large. (One specific 
example of this effect can be found in Allen & Koranda (1994); Koranda & Allen (1995). 
In those papers, we give a rigorous calculation of the effects of a stochastic gravity-wave 
background on the CMBR temperature distribution, in an inflationary model of the early 
universe.) 

If we combine the COBE observational limits on the different multipole moments 
{2 < I < 30) we arrive at the constraint on fig^^ given by 

,2 ^ -7 in-ii ( -"0 



f7gw/ifoo < 7 X 10"^' [-jj for Ho < f < SOHq. (4.5) 

Note that this limit does not apply to any gravitational waves, but only to those of 
cosmological origin, which were already present at the time of last scattering of the 
CMBR. This limit applies only over a narrow band of frequencies at the very lowest 
limits of Fig. |l|; it is far below any frequencies that will be accessible to investigation 
either by sp ace-b ased or earth-based experiments. 

In Section ^t] we will make a detailed examination of inflationary cosmological models, 
and investigate these limits on fJgw in more detail. 

4.2. Millisecond pulsar timing constraints on VLg^ 
For about a decade, Taylor and his research group (Kaspi et al. 1994)| have been mon- 



itoring the radio pulses arriving from a number of stable millisecond pulsars. Two of 
these pulsars are described in Table |l| These objects have turned out to be remarkably 
stable clocks; the regularity of their pulses places tight constraints on flg^ at frequencies 
of order 1/observationtime « lO^^Hz. 

The regularity of these pulsar's pulses can be described in terms of timing residuals. 
Roughly speaking, after fltting the arrival times of the pulses to a uniformly-spaced set 
(in the rest frame of the pulsars) the timing residuals are the remaining errors in the 
emission times of the pulses. As can be seen from Table the timing residuals are 
remarkably small; in dimensionless terms, 

^ « < lO-'\ (4.6) 

t 8 years 

It is easy to see that these extremely precise measurements put very tight restrictions on 
the amplitude of any gravitational waves with periods comparable to 8 years. 

To see the exact connection, imagine that a gravitational wave is passing by us, coming 
from a direction transverse to our line-of-sight to the pulsar, and assume its wavelength 
is 8 light years. Assume in addition that the pulsar is an ideal clock, emitting pulses 
at perfectly regular intervals in its rest frame. Then, during the time when one of the 
maxima of the gravitational wave is passing by the earth, the pulsar pulses will arrive 
spaced slightly farther apart (you can think of this as gravitational redshifting due to 
the gravitational wave). Four years later, a minimum of the gravitational wave will be 
passing by the earth, and the pulsar pulses arrive at an earth-bound observer slightly 
more closely spaced together (which can be thought of as gravitational blueshifting) . In 
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fact this system is just like a single-arm gravitational wave detector, with an arm length 
L of 8 light years. The size of the characteristic gravity- wave strain must be less than of 
order 

K< AL/L = At/t <10-^^, (4.7) 

where L = ct = 8 years and At ft! fisec is an upper bound on the residual tim- 
ing error. Substituting / = (8 years)""'^ = 10~^ Hz in Eqn. (|2.9[) gives the relation 



he = 10^^°-\/fi(/ ~ lO^^Hz) between the characteristic strain and the fractional energy 
density in gravitational waves. Hence the extremely small size of the fractional timing 
errors restricts he to be extremely small, which in turn gives a bound on f2 of 

n{f = lO^^Hz) < 10"**. (4.8) 

This bound provides important constraints on cosmic string networks, as discussed in 



more detail in Section 5.2. However it is not directly relevant for LIGO and other ground- 
based detectors, because the bound is at a frequency 10 orders of magnitude lower than 
their band of observation frequencies. 

The sensitivity of these pulsar timing experiments is truly remarkable! For example, 
the authors point out that in their 7-year data set for pulsar PSR B1855-I-09, the pulsar 
rotated precisely 40,879,349,533 times on its axis! They also conclude that remaining 
timing residuals are not due to gravitational waves: "We do not consider it likely that 
the observed trends in the timing residuals are in fact caused by gravity waves." 

I should note that a re-analysis of this pulsar timing data has recently appeared 



(Thorsett & Dewey 1996), which gives even more stringent limits (at 90% confidence) 
of r^gw < 4.8 X lO"'^ at 4.4 x 10~^ Hz. However this analysis has been criticized in 



other work |(McHugh et al. 1996)| which uses a different method (Bayesian rather than 



frequentist Ncyman-Pearson tests) to analyze the data, so it is still not entirely certain 
that the more stringent limit is reliable. 

4.3. Nucleosynthesis limits on fig^ 
The final bound on Q is obtained from the standard model of big-bang nucleosynthe- 



sis (often referred to as BBN (Kolb fc Turner 1990) ). This model provides remarkably 



accurate fits to the observed abundances of the light elements in the universe, tightly 
constraining a number of key cosmological parameters. One of the parameters which is 
constrained in this way is the expansion rate of the universe at the time of nucleosynthe- 
sis. This in turn places a constraint on the energy-density of the universe at that time, 
which in turn constrains the energy-density in a cosmological background of gravitational 
radiation. 

The BBN model's constraint can be expressed in terms of the number of massless 
neutrinos, which can be present and in thermal equilibrium at the time of nucleosynthesis. 



One typically obtains a bound of the form (Kolb & Turner 1990) 



Wn f \r) < |(iV,. - 3) / Prad \ , , 

/ / ^ ^gw at nucleosynthesis — ^ ^ 7io 7\ / V / 

l-|-3Xg+2Xg \ Pcrit / at nucleosynthesis 

where the quantities appearing on each side are the fractional energy densities in gravi- 
tational radiation and in relativistic matter at the time of nucleosynthesis. The integral 
includes only frequencies / greater than the Hubble expansion rate at the time of nu- 
cleosynthesis, in other words wavelengths smaller than the Hubble radius at that time. 
The BBN model constraints can be expressed in terms of the number of massless neu- 
trinos permitted; this is Nj, < 3.4. To give a rough idea of the limits on the present-day 
spectrum of Jlgw, we can convert this limit at the time of nucleosynthesis to a limit 
today. Since the energy-density of radiation (relative to the dominant energy density in 
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the universe) redshifts oc {1 + Z)~^ after the time of equal matter and radiation energy 
densities, this inserts a factor of (1 + Zeq)~^ into the limits. One obtains a bound on 
flgvf today of 



where Z^q « 6000 is the redshift at the time of equal matter and radiation energy 
densities. The range of integration now includes only wavelengths shorter than about 10 
pc. This bound restricts the spectrum of gravitational radiation Q^y, over a broad range 
of frequencies, but is not a very restrictive constraint. 



5. Cosmological sources of a stochastic background 

In this section, we will examine in detail three different physical processes which might 
have taken place in the early universe, which give rise to different spectra of flg„. In 
addition to their general interest, these three examples serve to motivate a number of 
interesting observations and comments. 

The first model we examine is the standard "inflationary" scenario for the early uni- 
verse (Kolb & Turner 1990; Lindc 1990). In this model, as the universe cooled, it passed 
through a phase in which the energy-density of the universe was dominated by vacuum 
energy, and the scale factor increased exponentially-rapidly. You will see that in its 
simplest version, this scenario is tightly constrained by the CMBR observational data 
(Allen & Koranda 1994; Koranda & Allen 1995). When this constraint is taken into 
account, the stochastic gravity wave background predicted by the simplest inflationary 
cosmological models are far too weak to be observable by either the initial or advanced 
LIGO detectors. 

The second model we examine is the cosmic string scenario (Vilenkin 1985; Vilenkin & 
Shellard 1994). In this model, as the universe cooled, long stringlike defects were formed 
at a phase transition. These "cosmic strings" form a network which self-intersects and 
chops off small loops of string. These small loops oscillate relativistically, and as they do 
so they emit a characteristic spectrum of gravitational waves. In contrast with the previ- 
ous example, for reasonable values of the parameters, the cosmic string scenario predicts 
a stochastic background which is large; enough to be observable; with the advanc;ed LIGO 
detector. In common with the previous example, the string network is described by a 
"scaling" solution, with the consequence that the spectrum of ilgw is flat over a wide 
range of frequencies. 

The third and final example in this section is also motivated by a phase transition in the 
early universe. We consider the case in which the phase transition is strongly flrst-order, 
and creates vacuum bubbles, wliic:li are bubbles of the new (low energy density) phase 
expanding in the old (high energy density) phase (Kamionkowski et al. 1994; Kosowsky 
et al. 1992; Kosowsky & Turner 1993). Within a short time, these bubble walls are 
moving relativistically, and much of the latent energy released by the phase transition 
is transformed into kinetic energy of the bubble walls. The highly-relativistic collisions 
of these bubble walls are a copious source of gravitational radiation, strong enough to 
be observable (in some cases) with the initial generation of LIGO detectors. Unlike the 
previous two examples, this phase transition/gravitational radiation production process 
is a "one-time" event, characterized by a particular cosmological time and a particular 
frequency today. Hence, in contrast with the previous two examples, you will see that 
the spectrum rigw(/) is strongly peaked at a frequency characteristic of this time. 

These three examples are intended primarily to be illustrative. They are by no means 




(4.10) 
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Initial (t=t^_ ) Final (t=t^ ) 




Figure 9. The pendulum consists of a mass hanging from a string. Initially the string is 
short, and the oscillation frequency is large. Between initial and final times ti and t/, the string 
is reeled out, increasing its length and decreasing the oscillation frequency of the mass. The 
quantum mechanical behavior of this system is more interesting than its classical behavior, and 
provides a simple model for the mechanism that generates large classical perturbations from 
zero-point fluctuations in inflationary cosmological models. 

comprehensive. The fact of the matter is that we don't know what happened in the 
early universe - the subject is extremely speculative, and will remain so until we have 
data (such as the spectrum of f2gw) that carries specific information about very early 
times. Nevertheless, these examples are quite useful, because they illustrate some of the 
mechanisms that might produce a cosmological stochastic background. They certainly 
lead mc to suspect that such a background is present, and to hope that it might be 
observable. 

5.1. Inflationary cosmological models 

Inflationary models of the early universe were studied extensively in the early- and mid- 
1980's (for a review, see Linde (1990)). Generally speaking, these are a class of cos- 
mological models in which the universe undergoes a phase of very rapid (power-law or 
exponential in time) expansion at early times. Inflationary models have a number of 
nice properties. They are very simple, and highly predictive. They provide "natural" 
solutions to the horizon, flatness, and monopole problems. They are also in good agree- 
ment with the COBE observations of the spectrum of temperature perturbations in the 
CMBR. (Note, however, that many other types of early-universe models are also in good 
agreement with these observations!) 

One of the remarkable facts about inflationary models is that they contain a beautiful 
mechanism which creates perturbations in all fields. These provide, within the context 
of inflationary models, a natural mechanism to create the density perturbations which 
evolved to form galaxies and clusters of galaxies today. More important for our purposes 
is that this mechanism also gives rise to an extremely distinctive spectrum of stochastic 
gravitational radiation. 

In this lecture, I will show you how these perturbations in inflationary cosmologies arise 
from the most basic quantum mechanical effect: the uncertainty principle. In effect, the 
spectrum of gravitational radiation that might be observed today is nothing other than 
adiabatically-amplified zero-point fluctuations. 

Let me begin with a simple example. Suppose that we have a harmonic oscillator like 
the one shown in Fig. 0. It is a traditional pendulum, consisting of a mass on a string. 
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In this case however the string is wound around a reel, and by unwinding the reel one 
can increase the length of the string and decrease the oscillation frequency ui = \fgjl 
of the system. Here / is the length of the string and g is the acceleration of gravity at 
the surface of the earth. We will denote the initial frequency by uji = y/g/h and the 
final frequency by = \/ sl^S- Let us suppose that the mass is initially at rest in its 
equilibrium position and consider in turn its classical and its quantum behavior. 

The classical behavior of this system is not very exciting. Since the mass is at rest in 
its minimum energy configuration (hanging straight down) it simply remains there. After 
reeling out the string, and lowering the frequency of the oscillator, the mass remains at 
rest, hanging straight down. 

The quantum mechanical behavior of this system is far more interesting. It has a long 
history (neatly encapsulated as the "adiabatic theorem" ) dating back to the invention of 
quantum mechanics in the 1920's. Because of the uncertainty principle, it is impossible 
for the oscillator to be truly at rest; in the minimum energy state it has initial energy 
Ei = ^huji, and executes zero-point motion about its classical rest position. In this initial 
state, the initial number of quanta Ni — 0. Quantum mechanically, this system has two 
possible types of behavior, depending upon how the string is reeled out. To distinguish 
these two cases, let T = tf —ti denote the time over which the string was steadily reeled 
out. 

The case where T is quite long compared to both the initial and final oscillation periods, 
T >> 27r/(jj/, is the slow or adiabatic case. In this case, the final energy of the oscillator 
is Ef = jTiojf. The difference between this energy and the initial energy represents work 
done to the reel that lowered the string. In addition to the "classical" work mgh done 
to the friction mechanism of the reel, a bit of additional "quantum" work Ei — Ef is also 
done. In this case, no quanta are created, TV/ = and the mass is still in its ground state 
after being lowered down. 

The truly fascinating case is the one in which T is quite short compared to the oscil- 
lation periods, T << 27r/wi. (Strictly speaking, we need l/l > cj at every moment.) In 
this case, the behavior of this quantum mechanical oscillator is quite different than the 
behavior of the corresponding classical oscillator. In this I will shortly show you, 

the final energy is half the initial energy: Ef — ^Ei = jhiOi. In this case, we create 
quanta: Nf = It is in this way that the inflationary universe models create large 

macroscopic fluctuations today, starting only with the minimum fluctuations expected 
from the uncertainty principle. In order to solve the horizon and flatness problems, in- 
flationary models require the universe to expand by a factor greater than w 10^^. This 
redshifts frequencies by the same amount, so that the ratio ^ > 10^^ and one ends up 
with large occupation numbers at the present. 

Before going into the details of inflation, let me show you a simple calculation to 
justify and explain the claims made in the previous paragraph. Shown in Fig. |l^ are the 
initial and final potentials of the harmonic oscillator pictured previously in Fig. ^ Let 
ai = mwi/h where m is the mass of the ball. The initial wavefunction corresponding to 
the ground state is 

^[x) = (5.1) 

Imagine now that at time t = this potential is suddenly changed to a final potential 
which is another harmonic oscillator potential but with much lower frequency wy, as 
shown in Fig. |l^. Let a/ = mtOf /%. Then we can express the wavefunction at time t = 
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Figure 10. A particle moving in the initial harmonic oscillator potential is in the ground state, 
and has a wavefunction shown by ip{x). The frequency of the oscillator is suddenly decreased, by 
external means, and the potential becomes that shown as the final potential. The wavefunction 
no longer corresponds to the ground state, but to a highly-excited many-particle state. 

in terms of the energy eigenstates of the new harmonic oscillator potential: 

oo 

ijjix) = ^ c„V„(a;) (5.2) 

where ipn{x) are the energy eigenfunctions of the final (lower frequency) potential: 

Mx) = NnHni^x)e-^''\ (5.3) 

Here Nn is a normalization constant, and the Hn are Hermite polynomials. In the usual 
way, we can work out the values of the expansion coefficients c„ : we multiply both sides 



of ( |5.2D by ^to(x) and integrate over all x, using the orthogonality of the eigenfunctions 



to obtain a formula for c„ 

Cn — for n odd, and (5-4) 



/ 2(2n)! 1 


( UJf -UJi \ 


/ ai + af n\ 





y ai + af nl \2{ujf + uji) J 
It is now straightforward to calculate the expectation value of the energy: 

Ef = f2{2n + l/2)hc.j\c,n\' = nu;f Q + ^^^J'^' ] . (5.6) 

n— ^i^f J 

It is easy to see that in the limit when the final frequency is much lower than the initial 
frequency, w/ << w^, one obtains Ef k, ^huji — hiOf {\^^ ■ This allows us to identify 

the number of created quanta as TV/ = (j^j • 

This simple mechanical model shows exactly how inflationary models generate pertur- 
bations in a "natural" way. Every mode of every quantum field can be thought of as 
a simple harmonic oscillator, and when the universe expands exponentially rapidly in 
size, the frequencies of these modes drop extremely rapidly. As you will see shortly, in 
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Figure 11. The cosmological scale factor a{rj) as a function of conformal time rj, for a 
simple inflationary model. The universe is deSitter space until time r] — rji when it makes an 
instantaneous phase transition to a radiation-dominated epoch. This would be a good "crude" 
model of our own universe, if slightly to the past of the present time rjo the universe became 
matter-dominated and a increased somewhat more rapidly between that time and the present. 



Table 2. The cosmological scale factor a{rj) for a two-stage inflationary cosmology. 

Time Phase Equation of State Scale Factor a{rj) Hubble H{ri) 

T] < rji DeSitter —P ^ p ^ const. rifrig^{2rii~r])~^ Hds = ci]o/Vi 

f] > rji Radiation P — p/3 v/vo cr)o/rf 



inflationary models, the universe expands rapidly in size, and because this expansion is 
very quick, it is not adiabaticl In other words, it corresponds to our second case above, 
where quanta are created. Even though this process is strictly quantum mechanical, and 
the energy created is proportional to h, the occupation numbers today are large, and 
these perturbations appear today to be large macroscopic fluctuations. In the case of 
the gravitational field itself, this process creates a characteristic spectrum of stochastic 
gravitational waves. Let me now sketch the calculation for a simple inflationary model. 
The details of this calculation can be found in Allen (1988). 

We assume that the universe is described by a simple "two-stage" cosmology as shown 
in Fig. |ll|. The metric of spacetime is given by 

ds^ = a^{r]){-drf + dx^ + hab{v, x)dx°'dx^), (5.7) 

where the scale factor a(ry), the Hubble expansion rate H{ri) and the equation of state 
are given in Table ^. (Note: ?7 is a conformal time coordinate with units of length.) Note 
that although the expansion does not appear to be very rapid during the deSitter phase 
(j] < r]i) this is because the scale factor 0(77) has been expressed in terms of a conformal 
time coordinate. When expressed in terms of physical or comoving time defined by 
cdt = a{r])dr] the scale factor is oc e^'*=* during the deSitter phase, and oc t^^"^ during 
the radiation-dominated phase. In order to solve the horizon and flatness problems. 
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enough expansion is needed to ensure that a{rjo) / a{rii) > 10^^ for GUT-scale inflation. 
So typically one has rjQ > IO^^t^i. 

The gravitational waves in this simple two-stage cosmology are described by metric 
perturbations of the form (Allen 1988) 



hab = eab{k)(f>{T]) exp(ifc • x), 



(5.8) 



where fc is a constant (in time rj) wavevector. The polarization tensors Cab are given 
in equations ( A 7 , A 8|) of Appen dix |A|. The a mplitude (j) satisfies a wave equation (the 
massless Klein-Gordon equation (Allen 1988)) 



^2 da d ^ 
drf a dt] drj 



= 0. 



(5.9) 



The solutions to this wave equation may be easily found in terms of elementary functions 
(which are simple cases of half-integer Bessel or Hankel functions) in both the deSitter 
and radiation-dominated epochs: 



For T] < rji (f){r]) — 
For rj > rji 4>(ri) — 



,-'ifc(»j-»7i) 



ae 



-1/0(77—771) 



+ pe 



ikiri—rii) 



(5.10) 
(5.11) 



where u — ck/a is the angular frequency of the wave, a and [3 are time- independent 
constants. Note that because k — \k\ is constant, the angular frequency w is a function 
of time. Since this differential equation is second order, one may obtain the values of the 
constants a and /3 by demanding that both (/> and d(f)/drj are continuous at the boundary 
rj = rji between the deSitter and radiation epochs of expansion. Imposing this restriction 
on the solution leads to 



1 

HnH, 



(5.12) 
(5.13) 



In these expressions, we use oj — ck/a^rja) to denote the angular frequency that would 
be observed today, and Hq = c/rjQ for the Hubble expansion rate that would be observed 
today. In effect, this problem is identical to tunneling problems done in elementary quan- 
tum mechanics, and we are simply calculating transmission/reflection coefficients. These 
types of calculations are also frequently referred to as Bogoliubov coefficient methods in 



the literature (Birrell fc Davies 1982) 



In an inflationary model, the long period of inflation damps out any classical or macro- 
scopic perturbations, leaving behind only the minimum allowed level of fluctuation, re- 
quired by the uncertainty principle. We have chosen our wave-function (5.10) to corre- 
spond precisely to this deSitter vacuum state. Provided that the period of inflation was 
long enough, the observable properties of the universe today should be indistinguishable 
from those of a universe which started in this vacuum state. We could work out the 
correct overall normalization of the wavefunction, but it is not necessary. 

Today, in the radiation-dominated phase of the universe, the eigenmodes describing 
particles are the coefficients of a, and those describing antiparticles are the coefficients 
of (3. One may easily show that the number of created particles of angular frequency u) 
today is given by (Allen 1988; Birrell & Davies 1982) 

= Number of created particles of freq uj = |/3ajp = — - — r^- (5-14) 
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Along the lines of the calculation given in Section 2.1 and in particular in analogy with 
(2.1), we may write an expression for the stochastic gravitational- wave energy-density 
contained in the (angular) frequency interval {(-u,U! + dio) as 

This may be re- written in terms of the present-day energy-density of the universe, and the 
(constant) energy-density during the deSitter phase, in the following way. The Hubble 
expansion rates are given in terms of these energy-densities by Hq = ^^^-f" and H^^ = 
^^^J'^' ■ The spectrum of stochastic gravitational waves is thus 

Pc df 9 & 9 ppianck 

where we have introduced the Planck energy density 

ppianck = ^ = 10"^ ergs/cm^ (5.17) 

You will notice that llgw is proportional to h, because this background is nothing other 
than (parametrically-amplified) zero-point energy! 

Several comments are now appropriate. First, our calculation of figw is for a simplified 
cosmological model that does not include a matter-dominated phase after the radiation- 
dominated epoch of expansion. If this matter-dominated phase is included, the result 



obtained is [(Allen 1988) 



^^gw(/) = ^^^(l + ^cq)-\ (5.18) 

y PPlanck 



for those waves which, at the time the universe became matter-dominated, had a wave- 
length shorter than the Hubble length at that time. Today this corresponds to frequencies 
/ > (1 -|- Zeq^^^Hf), where Z^q is the redshift of the universe when the matter and radi- 
ation energy densities were equal. At lower frequencies, the spectrum of figw(/) oc /^^. 
Second, note that our result (which is frequency-dependent) does not truly apply "from 
DC to light". For waves which are long compared to the Hubble length (10^^ cm today) 
corresponding to frequencies less than Hq today, the notion of energy density becomes 
non-sensical, because the wavelength becomes longer than the curvature length scale of 
the background space-time. In similar fashion, at high frequencies there is a maximum 
frequency above which f2gw drops rapidly to zero. In our calculation, we assume that the 
phase transition from the deSitter to the radiation-dominated epoch is instantaneous. 
However this process occurs over some time scale At, and above a frequency 

/max=-7-^-r- 5.19 

a{tQ) At 

which is just the redshifted characteristic rate of the transition, fig^ drops rapidly to 
zero. These low- and high-frequency cutoffs to the spectrum guarantee that the total 
energy-density in gravitational waves is finite rather than infinite. 

A typical spectrum of gravitational waves resulting from inflation is shown in Fig. p^ . 
For GUT energy-scale inflation, the ratio 

GUT - scale Inflation : — w (10^^ GeV/lO^^ GeV)^ w 10"^^ (5.20) 

PPlanck 

(Also shown for comparison is the naive 0.9 K blackbody, discussed in Section |2.l| .) 
We have also chosen /max ~ 10^° Hz, a reasonable value for GUT-scale inflation. You 
can see right away that in the inflationary case, the constraint coming from the COBE 
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Figure 12. The spectrum of stochastic gravitational waves in inflationary models is flat 
over a wide range of frequencies, and is shown as the solid curve. The horizontal axis is logj^g 
of frequency, in Hz. The vertical axis is logjpfJgw. The inflationary spectrum rises rapidly 
at low frequencies (wave which re-entered the Hubble sphere after the universe became mat- 
ter-dominated) and faUs off above the (appropriately redshifted) frequency scale /max associated 
with the fastest characteristic time of the phase transition at the end of inflation. The amplitude 
of the flat region depends only on the energy- density during the inflationary stage; we have cho- 
sen the largest amplitude consistent with the COBE constraint discussed earlier: f2g„ < 7x10"^^^ 



at 10 Hz. This means that at LIGO and LISA frequencies, f2gw < 8 x 10 
comparison is the spectrum of a 0.9 K blackbody, as discussed earlier. 



Shown for 



observations discussed in Section |4.l| puts very severe restrictions on the spectrum (Allen 
& Koranda 1994; Koranda & Allen 1995). In Fig. |l^ we have chosen the amplitude of 



the spectrum (determined by the ratio 



PPla 



-) to be as large as possible, consistent 



with this COBE constraint. Nevertheless, because the spectrum falls off cx /^^ at low 
frequencies, this means that today, at LIGO and LISA frequencies (indicated by the 
arrows), ilgw < 8 x 10^^"*. This is certainly too small to be observable with either initial 



or advanced LIGO, 
LISA. 



and is at the limits of what might be potentially observable with 



Recent work (Turner 1996) has examined in more detail the spectrum of stochastic 
gravitational radiation produced in slow-roll inflationary models. In these models, the 
spectrum is not exactly flat at high frequencies, but has a (negative) logarithmic slope 
riT- There is a connection between this slope and the amplitude of the spectrum at the 
lowest frequencies; if the amplitude is large then nx is quite negative, and if ut is near 
zero then the amplitude is small. Consequently Turner found that in slow-roll inflation 
the upper bounds on Qg^ at LIGO and LISA frequencies are about an order of magnitude 
lower than the ones given here. 

There are other types of "inflationary" cosmology based on superstring models of 
the fundamental interactions, which have recently received a great deal of attention 
(Brustein et al. 1995). These models typically contain a "deflationary" epoch followed 
by a stage of rapid expansion, and produce a spectrum in which r2gw(/) oc below 
some maximum cutoff frequency. This spectrum evades the CMBR and millisecond 
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pulsar timing constraints, yet might still produce enough stochastic background to be 



observable by LIGO (Allen fc Brustcin 1996) 



Before moving on to our next example (cosmic strings) of how cosmological processes 
can give rise to a stochastic background, a few final comments are in order. The first 
is that the long flat region of the spectrum rig„ reflects the scaling behavior of the 
cosmological scale factor a{t) (x exp{H^st) during the deSitter phase of expansion. We 
will see identical self-similar behavior in our next example. Notice that the highest 
frequency parts of the spectrum (such as the falloff at lO^*' Hz in our example) reflect the 
behavior of the universe at early times (in this example, the nature of the phase transition 
that ended inflation at « 10~^^ seconds after the big bang). The low frequency parts of 
the spectrum reflect the fairly recent history of the universe. For example, for inflationary 
cosmology, the "kink" in the spectrum at 10^^^ Hz reflects the change in the expansion 
law from radiation- to matter-domination at Zcq ~ 6000. Grishchuk has shown that 
in general, this type of mechanism (parametric amplification of zero-point fluctuations) 
produces a spectrum of the sort we have illustrated. However if the universe does not 
expand exponentially (constant H{t)), but "speeds up" and "slows down" (i.e. H{t) 
increases and decreases) then this is reflected in the spectrum of Qg^j which is no longer 
flat, but instead has characteristic increases and decreases, reflecting the changing values 
of H{t) during the evolution of the universe. 

5.2. Cosmic string cosmology 

In this lecture, I am going to show you a calculation of the stochastic gravitational 
wave background produced by a network of cosmic strings (Vilenkin 1985; Vilenkin & 
Shellard 1994). These are objects which might have formed during a phase transition, as 
the universe cooled. We do not know if our universe contains cosmic strings or not (this 
depends upon the precise nature of the phase transition) so these calculations must be 
viewed as illustrative rather than realistic. 

Cosmic strings are one-dimensional (in space) string-like objects. They are topologi- 
cally-stable scalar/gauge field configurations, whose dominant decay mechanism is the 
emission of gravitational waves. These string-like objects are analogous to the vortex 
lines which form in superfluid helium, or to topological defects which can form during 
phase transitions of liquid crystals. Here, our interest is in cosmic strings which might 
have formed during a phase transition as the universe cooled, and a fundamental C/(l) 
local gauge symmetry was broken. These are the simplest strings; they are characterized 
by a single dimensional scale: their mass-per-unit-length fi. A phase transition which 
formed strings at the Grand Unified Theory (GUT) energy scale of about 10^^ GeV 
would result in strings with ^ « 10^^ gms/cm. 

These kinds of cosmic string have several remarkable properties. First, they are formed 
without any "ends" so the strings are always in the form of loops. Sometimes the loops 
are small (compared to the Hubble length) in which case we call them "loops" . Sometimes 
the loops are large (compared to the Hubble length) in which case we call them "infinite 
strings" . In neither case are there any unjoined "ends" . Note that in assigning a "length" 
to the loop, we are using an invariant measure of the length; the length is defined as the 
total energy of the loop divided by its mass-per-unit-length fi. A second remarkably 
property of these loops is that they have a tension equal to their mass-per-unit-length 
fi. This is in contrast to the kind of strings we use to wrap packages, which have small 
tension compared to their mass-per-unit-length /i. Because of this large tension, the 
cosmic strings oscillate relativistically under their own tension. Thus, a circular loop 
of cosmic string, initially placed at rest, will collapse under its own tension in a time 
comparable to the time it takes a light ray to move a distance equal to the radius of 
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Figure 13. A snapshot of a universe containing a cosmic string network. The spatial volume 
shown here is a cube with side length 2c/H where H is the Hubble expansion rate. The 
string network has a remarkable scaling property: at any time the appearance of the network is 
(statistically) like the one shown here, provided that all length scales are expressed in terms of 
the Hubble length at that time. The network consists of many tiny loops, which here appear as 
points, and a number of so-called "infinite" strings. Although these infinite strings appear to 
have ends, this is just an artifact of having "cut out" a Hubble volume cube for this snapshot. 



the loop. In general, a non-circular loop will oscillate quasi-periodically, with a period 
equal to (half) the time it takes light to travel a distance equal to the length of the 
string making up the loop. A third remarkable property of these cosmic string loops is 
that they are stable against all types of decay apart from the emission of gravitational 
radiation. So typically a loop of cosmic string, under the effects of its own tension, will 
oscillate almost periodically, gradually emitting gravitational radiation, and shrinking in 
size. 

One might worry that the configuration of these cosmic strings is very complicated, 
depending precisely upon how the mass-energy of the universe was distributed when 
they were formed, and that for this reason, it would be terribly difficult to make any 
precise predictions about the effects they would have on the universe. Fortunately, about 
fifteen years of extensive study of networks of cosmic string has revealed that they have 
a remarkable "scaling" property (Vilenkin 1985; Vilenkin & Shellard 1994). This can 
be summarized by Fig. |l^, which shows a snapshot (at one instant of cosmic time) of 
a universe containing cosmic strings. The crucial point is this one: at any instant in 
time, the appearance of the string network is (statistically) like that shown in the figure 
Fig. As the universe expands, the string network evolves; however at any time, 
each Hubble sized volume contains, on the average, a certain number (say A=52) of 
infinite strings passing through it, and a very large number of small loops. The small 
loops oscillate, emitting gravitational radiation, and as the universe evolves, more small 
loops are chopped off the infinite strings, replacing those loops which are disappearing 
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Figure 14. A typical spectrum of gravitational radiation produced by a cosmic string network 
is shown as the solid curve. The dotted curves show the contributions of loops which "died" 
while the universe was radiation-dominated, and the contributions of loops which "died" more 
recently, when the universe was matter-dominated. The part of the spectrum which lies in the 
LIGO band comes from the radiation-dominated epoch. Also shown, for naive comparison, is 
the spectrum of a 0.9 Kelvin blackbody. The frequency /obs is where the millisecond pulsar 
timing observations constrain the spectrum. 



after converting their energy into gravitational waves. All of the properties of the network 
"scale" with the Hubble length. So even though it is impossible to predict the appearance 
of a string network in our own universe, in a statistical sense, we can describe the average 
properties of the network at any time. 

The spectrum of gravitational radiation produced by a cosmic string network has 
been carefully studied (Vilenkin 1985; Vilenkin & Shellard 1994; Caldwell & Allen 1992; 
Caldwell 1993). The spectrum of gravitational radiation produced by a "typical" cosmic 
string network is shown in Fig. You can see that although the spectrum is not 
flat, it is fairly flat in the range of frequencies which LIGO and the other ground-based 
detectors will be sensitive to. This radiation was produced by loops which emitted their 
gravitational waves during the epoch when the universe was still radiation-dominated. I 
want to show you an estimate of the amplitude of the stochastic spectrum i7g„(/) for the 
range of frequencies / over which fig^ is roughly constant. To do this, I need to carry 
out counting arguments, to count the number of loops emitting gravitational radiation at 
any time. The number of loops that we have depends upon the volume that we consider, 
and the counting gets tricky if we don't make the right choice. For example, we could 
choose a volume which at any time was a single Hubble volume (a ball of radius equal to 
the Hubble length), or we could pick a constant co-moving volume, or a fixed physical 
volume. To avoid the need to make one of these choices, there is a simpler way to do the 
counting. We'll do our calculations for a universe which has a finite spatial volume, and 
whenever we count something (say loops), we'll simply count the total number of those 
things (i.e. loops) in the entire universe. 

The universe we will use is shown in Fig. HS. At any instant in time, its spatial 
topology is that of a three-dimensional torus T^— x x S^. The metric is ds'^ = 
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Figure 15. A spatial slice of the model universe, at one instant in time, is a three-dimensional 
torus. In the figure, we have shown a single long cosmic string, which winds around this 
universe 5 times in the a;— direction and 1 time in the y— direction. Because the figure shows 
a two-dimensional projection, we can't see how many times the string winds around in the 
2— direction. Note that the string has no ends. Also shown is a small sphere, which is the 
boundary of the Hubble volume at the given instant in time. At later times, this small sphere 
would be larger, and would include a larger part of the three-torus. The single long string shown 
in this figure is a so-called "infinite" string. Although its length is finite and it has no ends, it 
is not called a "loop" because it is longer than the Hubble length. 

—c^dt^ -t- a?{t)dSP'^ but the spatial coordinates are each periodically identified on the 
interval [0,L]. Thus each of the three spatial coordinates ranges over < x,y, z < L. 
The time coordinate t is physical or comoving time, and during the radiation-dominated 
epoch of cosmological expansion, a{t) — {t/toY^'^ where tg is the present time, related 
to the present Hubble expansion rate by to = {2Hq)~^. Our model universe has a 
finite proper spatial volume; at any instant in time this three volume is V{t) = L^a^{t). 
Provided that this volume is greater than the Hubble volume (in other words, t is less 
than some critical value) then the fact that our model universe is finite in spatial extent, 
rather than infinite, does not matter and observable quantities will have the same values 
that they would have in a truly infinite universe. So the toroidal topology of our model 
universe simply serves as a convenient calculational tool. 

Shown in Fig. |l^ is the spectrum of gravitational radiation produced by a cosmic string 
network described by a particular set of dimensionless parameters. This spectrum was 
actually calculated by a computer program (Caldwell & Allen 1992; Caldwell 1993); there 
is no way to write a useful analytic expression for this spectra. You will notice however 
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Figure 16. The spectrum of gravitational radiation produced by a cosmic string network. 
Around the LIGO frequency band / ~ 100 Hz, the spectrum is fairly flat. In these lectures, we 
derive a simple analytic estimate for the amplitude of Qg„ in this flat region. 



that the LIGO frequency band (around 100 Hz) ilgw is flat (a constant, independent of 
frequency). It is not difficult to derive a formula for this constant value. 

To analytically determine the amplitude of Qgw in the LIGO band, we will make 
extensive use of the scaling property of the string network. This property has been well 
established, both by computer simulation and by analytic arguments (Vilenkin 1985; 
Vilenkin & Shellard 1994). The scaling property may be neatly summarized: A single 
correlation length l{t) = Hubble length = 2ct characterizes ALL properties of the string 
network. This scaling property can be used to obtain qualitative relations. For example, 
the energy-density of long strings is given by 

Afilc'^ An 
= = (5-21) 

where pL is the mass-per-unit-length of the string, and ^ is a dimensionless constant 
(determined by numerical simulations to have a value A « 52). Roughly speaking, A 
is the number of long strings which pass through the Hubble volume shown in Fig. f^. 
(An example of such a long string passing through a Hubble volume is shown in Fig. |15| .) 
Another example of scaling concerns the sizes of loops of string which are "chopped off' 
the infinite string network. From scaling, the size of these loops at the time of formation 
must be 

Size of loop formed at ibirth = acibirth, (5.22) 
where a is a dimensionless constant (which in principle can be determined by numerical 
simulation) . 
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Because the gravitational radiation produced by a cosmic string network comes from 
the small oscillating loops of string, we need to determine the rate at which these loops 
are formed. This can be found from energy conservation. The decrease in the total 
energy of long strings in the universe must be compensated by a corresponding increase 
in the total energy of cosmic string loops: 

d{Vp^) + (^ate3)d7Vioops = 0. (5.23) 

(Note that the equation of state of the string loops is similar to that of dust; otherwise 
the equation of energy conservation would be somewhat more complicated.) Using the 
expressions above for and for V we obtain an equation for the rate at which loops of 
string are formed: 

^^l°°P'^ _ 1 ,-3/2 5/2 

where Tve used the explicit time dependence of a{t) . The loops which are being cut off the 
infinite strings oscillate relativistically, emitting gravitational radiation. A remarkable 
fact about these loops is that the energy-loss rate (power radiated) P does not depend 
upon the size of the loops, but only on their shape (Vilenkin 1985; Vilenkin & Shellard 
1994). For an "average" loop 

P = "fGfi^c, (5.25) 

where 7 is a dimensionless constant, which has been determined numerically to have a 
value 7 « 50. Because the loops radiate power at a constant rate, independent of their 
size, the energy and mass of a loop decrease linearly with time, until all the energy has 
been radiated away. At time t, a loop formed at time tbirth has length 

(Size of loop formed at ^birth, at time t) = actbirth — lG^{t — tbirth)/c. (5.26) 

The loop disappears when this length reaches zero, at a time 

ideath = 1 ^ I ^birth = /3ibirth, (5-27) 



where we have defined a constant (3 — 1 + > 1 which is the (constant) ratio of 
death-time to birth-time for cosmic string loops. We are now in a position to work out 
an analytic approximation to Qg^{f) for cosmic string networks, which is valid in the 
LIGO band. 

Consider the total energy (as seen today) emitted in gravitational waves by all oscil- 
lating cosmic string loops which were formed (1) after the cosmic string network came 
into existence at time ^formation and (2) which are no longer present, and thus were born 
before time to / (3. This total energy is 

Eg^ = / dt'—-^ / dt"jGfi^c^---j (5.28) 

-I- (Ignore Loops Still Present Today). (5.29) 

The different factors appearing in this formula are (1) the number of loops diVioop born 
in the time interval dt' , (2) the energy dE = dfjCfi^c emitted by those loops in the time 
interval dt", and (3) the redshift factor a{t") / a{tQ) which reduces the energy as observed 
today. Substituting in expression (5.24) for the rate of loop formation and carrying out 
the integral over t" gives 

''""^ \l^G^lW''t'''' {P"' - l)l (5.30) 



dt'- 



dt' 
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Because of the scaling nature of the string network, and because frequencies redshift 
as the logarithmic frequency interval over which the gravitational radiation is 

distributed is precisely 



\ J low / \ '^formation / 



(5.32) 



The In /? term can be neglected in comparison with the other In term, because typically 
the magnitude of 100 < (3 < 10^ is much smaller than the magnitude of to/iformation ~ 
10^° for GUT scale strings. If we now use the relationship between the Hubble expansion 
rate and the critical energy density, 

„2 _ 1 _ 87rG/3critical /(. 

^0 - 4f - 3^^ ^^■^^> 
we can immediately calculate the spectrum of gravitational waves: 

f^gw = ^— ^ (5.34) 

Pcritical din J 

^critical ln(/high//iow) 

Now this calculation is for a radiation-dominated universe. If we include the phase of 
matter domination following it, the amplitude of the spectrum is decreased by a factor 
of (1 + Zcq) to give our final result for the amplitude of the stochastic gravity wave 
background of a cosmic string network, in the LIGO band: 

2 



_ 167r ^7 fG/j, 



K-l)(l + ^e.)-. (5.37) 

This quantity can be large enough to be observable with the LIGO detectors, for values 
of the dimensionless parameters characterizing tlw^ cosmic strings. 

The values of these dimensionless parameters are not completely known. Reasonable 
values are A = 52 for the number of long strings/horizon volume, 7 » 50 (for the 
radiation rate from a typical loop) and mass-per-unit-length Gfi/c'^ ~ 10^® for GUT 
phase transitions. However the value of a (the size of a loop at formation) is still 
unknown. The high-resolution numerical simulations have established that a < 10~^ 
but this is only an upper bound. Indeed it is thought that a could perhaps be as small 
as ^Gn/c^ 10~^, which is the scale at which gravitational back reaction cuts off the 
small scale structure and fluctuations on the long strings (Vilenkin 1985; Vilcnkin & 
Shellard 1994). This uncertainty in the value of a leads to a wide uncertainty in /3; the 
allowed range is /3 e [1,200]. Because this parameter governs the lifetime of a loop, it 
has a large effect on the spectrum flg^. 

There is increasing evidence that the scale of loop formation (i.e. the value of a) is 
indeed set by gravitational back reaction, and is thus very small. One can work out the 
limit a ^ of figw, obtaining 

I4f 



^^gw^^^H^ (l + ^eq)-'. (5.38) 



This may well turn out to be the correct expression in the LIGO frequency range. How- 
ever the allowable range of is tightly constrained by the millisecond pulsar timing 
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Figure 17. A first-order phase transition can produce rapidly-expanding bubbles containing 
the new (low temperature) phase within the old (high-temperature) phase. The diflerence in 
energy (latent heat) is transformed into kinetic energy of the bubble walls. When the bubbles 
collide, a fraction of this energy is emitted as gravitational radiation. 



observations discussed in Section 4.2, and the increasingly strong constraints of these 
observations may rule out GUT scale cosmic strings in the near future (Caldwell & Allen 

1992; Caldwell 1993). 

Recent work [(Martin fc Vilenkin 1996) has shown that "hybrid" topological defects 
could evade the increasingly strong constraints from pulsar timing observations. For 
example a cosmic string network might form at early times, as described in this section, 
but then at later times, domain walls form (having the strings as their boundaries) 
and "pull together" the strings, destroying the network. In such a scenario there is a 
high-frequency stochastic background (produced by the strings at early times) but no 
low-frequency gravitational radiation. In this way, the scenario evades the millisecond 
pulsar timing constraints on figw at lO^'* Hz. The prospects for detecting hybrid defects 
of this type with initial LIGO do not appear good, however advanced LIGO should have 
more than enough sensitivity. 



5.3. Bubbles from first-order phase transitions 

My final example of a specific mechanism that can produce a stochastic background of 
gravitational radiation is rather different than the two previous examples I have examined 
(inflationary cosmology and cosmic strings). In these two previous examples there was 
a scaling law or behavior which resulted in a spectrum f2gw(/) that was flat, at least for 
waves which entered the Hubble sphere before the universe became matter-dominated. 
In this final example, I will briefly review the effects of a different type of mechanism. 
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Figure 18. Typical spectra ngw(/) produced from the collision of bubbles which result from a 
first-order phase transition. This graph is reproduced from Fig. 7 of Kamionkowski et al. 1994. 
The horizontal axis is logjQ(27r///3) and the vertical axis is logj^g Slgw/3^(1 + a)^ H~^or^ . 
The spectrum peaks at a characteristic frequency /max, characteristic of the expansion rate and 
time at which the bubble collisions occurred. The spectra are shown for bubble wall velocities 
oiv = 0.2,0.4,0.6 and 1.0. 

In this mechanism, as the universe expands and cools, the matter within it undergoes a 
first-order phase transition at some definite time in the past (Kamionkowski et al. 1994; 
Kosowsky et al. 1992; Kosowsky & Turner 1993). This transition occurs when the tem- 
perature of the universe has dropped sufficiently below the characteristic energy of the 
phase transition: kT^ < -Btransition- Bubbles of the new (low-energy-density) phase form 
within the old (high-encrgy-density) phase, as shown in Fig. These bubbles expand 
rapidly after formation, converting the difference in energy density x bubble volume 
into kinetic energy of the bubble walls. Within a short time, these bubbles are moving 
relativistically, and within a Hubble expansion time, they collide. These collisions are 
highly relativistic, non-symmetric events, and produce copious amounts of gravitational 
radiation. Within a short time after the collisions begin, a significant fraction of the 
energy that was in the bubble walls is converted into gravitational radiation. 

In contrast to our previous two examples, this bubble collision process produces a 
spectrum which is strongly peaked at a particular frequency, characteristic of the time 
at which the phase transition and bubble collisions took place. This spectrum has been 
computed (Kamionkowski et al. 1994; Kosowsky et al. 1992; Kosowsky & Turner 1993) 
and is shown in Fig.|l^. The spectrum is peaked at a characteristic frequency 



The Hubble expansion rate at the time of the phase transition is denoted and /3 is a 
measure of the bubble nucleation rate. For typical phase transitions, Kosowsky, Turner 
and their collaborators (Kamionkowski et al. 1994; Kosowsky et al. 1992; Kosowsky & 
Turner 1993) find that f3/H^ « Aln^Mpi^nckC^ / kn) ~ 10^ The number of degrees of 
freedom is also of order 10^ in typical GUT models. For example, in the minimal 




(5.39) 
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standard model electro-weak phase transition one has the transition at 
102 GeV with /max « 4.1 X lO^^ Hz. 

The amplitude of the spectrum depends mostly upon the relative energy-density differ- 
ence between the interior and the exterior of a bubble. This difference in energy density 
is deposited directly in the bubble walls, determining the speed at which the walls crash 
together, and thus the amount of energy converted to gravitational radiation. The dif- 
ference is described by a parameter a = ^^1^™° ^ . Kosowsky, Turner and Watkins find 
that the amplitude of the spectrum at the peak is approximately 

n,M.,A 1.1 X io-„» (A)" (^)' (j^) (Ij)-'"^ (5,40) 

The parameter k is the fraction of vacuum energy that goes into kinetic energy of the 
fluid; it is an increasing function of a and typically lies in the range 10"^ < k < 1. The 
final parameter v is the propagation velocity of the bubble walls. For strongly first-order 
phase transitions, a — > oo and u — > 1. In this strongly first-order case, a phase transition 
at « 10^ GeV could be easily observed with the Advanced LIGO detectors. 

We should note that it is perhaps overly optimistic to assume that the phase transition 
is strongly first-order. At least one phase transition which is expected to occur, at an en- 
ergy scale of about 10^ GeV, is the electro- weak transition. In the standard SU(3) model 
this phase transition has 17gw(4.1 x 10"'^ Hz) k, 10"^^, which is many orders of magnitude 
too small to be observable by any proposed space-based detector (Kamionkowski et al. 
1994; Kosowsky et al. 1992; Kosowsky & Turner 1993). However it is perfectly possible 
that there as a strongly first-order phase transition in the earlier history of the universe, 
that might give rise to waves observable with the proposed LISA detector. 



6. Conclusion 

I am often asked, after a series of lectures like this, "Do you really think that a stochas- 
tic background is there?" There is only one possible honest answer ~ I do not know! The 
problem is that we have very little information about what really happened in the very 
early universe. And the scenarios that we have given as examples are just that, examples. 
Personally, I think it is unlikely that any of them will turn out to be true. However they 
do have illustrative value. They show that, under reasonable assumptions about what 
happened in the early universe, it is possible to produce a stochastic background which 
is easily detected with the instruments now under construction. The possibility of such 
detection has been anticipated for some time (Grishchuk 1976; 1988), and there are good 
hopes that they will be realized in the next decade. 

So while the answer to the question is unknown, that does not matter. In the next 
decade we will probably know the answer in a very direct way, with real observational 
data. And if we do detect a stochastic background, a new era in our understanding of 
the early universe will begin. 
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Appendix A. Calculational details 

The purpose of this section is to provide details of the calculations of the expectation 
values of the strains due to a stochastic background. The starting point is the plane 

wave expansion of the gravitational metric perturbations. In transverse traceless gauge, 
this can be written in the form of a plane wave expansion 

Kb{t,x)=Y, I d'kCA{k)sm{c\k\t-k-x + ^{k))ei^(k). (Al) 

A ■' 

Here A = +, x label the two polarization states, k is the wave-vector of the perturbation, 
Ca arc arbitrary real functions specifying the amplitude of the modes, $ is a real function 
taking values in the interval [0, tt] and e^^ are the polarization tensors, which depend only 
upon the direction of the wave-vector k. This decomposition is unique: every gravity 
wave can be represented exactly one way, in this form. 

Rather than representing the traveling wave by a phase and a real amplitude, it is 
more convenient to use a complex amplitude. An equivalent representation is 

hab{t, x)=^(^J^ df j^^ dh BAif, Cl) exp(27ri/(t - Q ■ x/c))e^,{n)^ (A 2) 

Here, n is a unit vector on the two-sphere; we have effectively written k = 2'KfCl. The 
function BA{f,^) is an arbitrary complex function of the two variables, defined on the 
range / > 0. By taking the real part of this expression explicitly, one can cast it in the 
form 

habit,x)=y^[ df j dClhA{f,n)eM'2TTifit-n-x/c:))e^t,(n), (A3) 

where now /ia(/; ^) is an arbitrary complex function satisfying the relation hA{—f, ^) = 
h\{f,fl). As before, this decomposition is miique. 

The polarization tensors appearing in these relations may be given explicitly. In stan- 



dard angular coordinates {9, (j)) on the sphere one has 

fi = cos sin + sin (/) sin 0y + cos 0£ (A 4) 

rh = sin (j)x — cos 0y (A 5) 

h = cos (p cos 6x + sin (p cos 9y — sin 6z (A 6) 

^abi^) = '^a.mb - UaUb (A 7) 

^abi^) = + 'f^a'mb (A 8) 



One can verify by inspection that m and fi arc a pair of orthogonal unit-length vectors 
in the plane perpendicular to Cl. It is simple to show that any rotation of the vectors rh 
and h within the plane that they define simply corresponds to a trivial re-definition of 
the complex wave amplitudes and h-^. 

We are now in a position to verify some of our key formula. We begin with the 
assumption that the stochastic background is stationary and Gaussian. This means that 
the ensemble average of the Fourier amplitudes is 

{h*A{f,h)hA'{f'A'))=5{f-f')5\h,(l')5AA'H{f) (A 9) 

for all / and /'. In these formula, 

5'^{Cl,Cl') = (5((/)- (/-')'5(cos6'- cos6'') (A 10) 

is the covariant Dirac J-function on the two-sphere, and the covariant volume element is 
dO, = smOdOdcj). The function H{f) is related to the spectrum r2gw(/). H{f) is a real. 
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non-negative function, satisfying H{f) — H{—f). To see the relationship, consider the 
energy density in gravitational waves. This is given (locally) by 

Pgw = 3^(/la6/i'^''), (All) 

where the overdot denotes a time derivative, and both tensors are evaluated at the same 
space-time point {t,x). Substituting the plane wave expansion into this formula yields 

/oc /• 
df dhATT^fH{f)ei,{h)et{h). (A 12) 

/I '°° •'^'^ 

Since Gab^°A — ^ ^^'^ J dd ^ An one has 

/OO POC 
df fH{f) - 1287r3 / df f^H{f). (A 13) 

-OO Jo 

Using the definition S7gw one obtains the relationship between the spectrum f2gw and 
H{f). For / > one has 



This formula is needed to obtain the next result. 

The second relation required is the expected value of the Fourier amplitudes at two 
different sites. This is now easily obtained. By definition 

{mMf')) 

= f dn f dCi'{h*Af,Ci)hA'{f',n')) 

AA' ■'S' 

X exp(27ri/f2 • xi/c-2mf'h' ■ X2I c)F^{h)Ff {(l') 

= V / dhH{f) exp(27rz/f^ • / c)F^{n)F^{n)6{f - /') 
A 

= Y^if)H{f)S{f f) = ^\f\-^n^^{\f\)^{\f\)5{f - /'). 

In the third line above, we have used the definition of the strain in the detector, given 
in ( 3.10| ), and in the fourth line, we have used the definition of the overlap reduction 
function (Eol). 



REFERENCES 
Abramovici, A., et al. 1992 Science 256, 325. 
Allen, B. 1988 Phys. Rev. D 37, 2078. 
Allen, B., Koranda, S. 1994 Phys. Rev. D 50, 3713. 

Allen, B., Brustein, R . 1996 Detecti ns relic gravitational radiation from string cosmology 
with LIGO, preprint 5r-qc/9609015 , submitted to Physical Review D. 

Allen, B., Ottewill, A. C. 1996. Detectio n of anisotropies in the gravitational-wave stochas- 
tic background, preprint gr-qc/9607068, submitted to Physical Review D. 

Apostolatos, T. a. 1995 Phys. Rev. D 52, 605. 

Apostolatos, T. a., Cutler, C. Sussman, G. J., Thorne, K. S. 1994 Phys. Rev. D 49, 
6274. 

Bennett, C. L., et al. 1992 Astrophys. J. 396, L7. 
Bennett, C. L., et al. 1994 Astrophys. J. 436, 423. 



44 B. Allen: Cosmological sources 

Bennett, C. L., et al. 1996 Four Year COBE DMR Cosmic Microwave Background Observa- 
tions: Maps and Basic Results, astro-ph/9601067, submitted to Astrophys. J. Lett.. 

BiRRELL, N. D., Davies, P. C. W. 1982 Quantum fields in curved space, Cambridge University 
Press, Cambridge. 

Blair, D., Lu, L. 1996 A cosmological background of gravitational waves produced by super- 
novae in the early universe, to appear in Mon. Not. R. Astron. Soc. 

Blanchet, L., Damour, T., Iyer, B. R., Will, C. M., Wiseman, A. G. 1995 Phys. Rev. 
Lett. 74, 3515. 

Brustein, R., Gasperini, M., Giovannini, M., Veneziano, G. 1995 Phys. Lett. B 361, 45. 

Caldwell, R. R. 1993 in Proceedings of the fifth Canadian General Relativity and Gravitation 
Conference, eds. R. McLenaghan and R. Mann, World Scientific. 

Caldwell, R. R., Allen, B. 1992 Phys. Rev. D 45, 3447. 

Caron, B., et al. 1995 in Gravitational Wave Experiments, proceedings of the Edoardo Amaldi 
Conference, World Scientific, 86. 

Christensen, N. 1992 Phys. Rev. D 46, 5250. 

Danzmann, K., et. al, 1995 in Gravitational Wave Experiments, proceedings of the Edoardo 
Amaldi Conference, World Scientific, 100. 

Dhurandhar, S. v., Sathyaprakash, B. S. 1994 Phys. Rev. D 49, 1707. 

Flanagan, E. 1993 Phys. Rev. D 48, 2389. Note that the second term on the r.h.s. of equation 
(b6) should read — lOji(Q) rather than — 2ji(q), and that the sliding delay function shown 
in Figure 2 on page 2394 is incorrect. 

Grishchuk, L. 1976 Sov. Phys. JETP Lett., 23, 293. 

Grishchuk, L. 1988 Usp. Fiz. Nauk., 156, 297 [Sov. Phys. Usp. 31, 940]. 

Hawking, S. W., Ellis, G. F. R. 1973 The Large Scale Structure of Spacetime, Cambridge 
University Press, Cambridge, 141. 

Hough, J., et al., 1995 in Gravitational Wave Experiments, proceedings of the Edoardo Amaldi 
Conference, World Scientific, 50. 

Kamionkowski, M., Kosowsky, a.. Turner, M. S. 1994 Phys. Rev. D 49, 2837. 

Kaspi, v., Taylor, J., Ryba, M. 1994 Astrophys. J. 428, 713. 

KOLB, E. W., Turner, M. 1990 The Early Universe, Frontiers in Physics, Addison Wesley. 
KORANDA, S., Allen, B. 1995 Phys. Rev. D 52, 1902. 
Kosowsky, A., Turner, M. S., Watkins, R. 1992 Phys. Rev. D 45, 4514. 
Kosowsky, A., Turner, M. S. 1993 Phys. Rev. D 47, 4372. 

LiNDE, A. D. 1990 Particle Physics and Inflationary Cosmology, Harwood, Chur, Switzerland. 

Martin, X., Vilenkin. A. 1996 Gravitational wave background from hybrid topological defects, 
preprint astro-ph/9606022, to appear in Physical Review Letters. 

McHuGH, M.P., Zalamansky, G., Vernotte, F.,Lantz, E. 1996, preprint, Besangon Ob- 
servatory. 

MiCHELSON, P. F. 1987 Mon. Not. Roy Astron. Soc. 227, 933. 
NiEBAUER, T.M., et al. 1993 Phys. Rev. D 47, 3106. 

Owen, B. 1995 Search templates for gravitational waves from inspiraling binaries: choice of 
template spacing, Caltech preprint, QC-9511032, submitted to Physical Review D. 

Sachs R. K., Wolfe A. M. 1967 Astrophys. J. 147, 73. 

SCHUTZ, B. F. 1991 in The Detection of Gravitational Waves, ed. D. G. Blair, Cambridge 
University Press, 406. 

Smoot, G. F., et al. 1992 Astrophys. J. 396, LI. 

Tagoshi, H., Sasaki, M. 1994 Prog. Theor. Phys. 92, 745. 

Thorne, K. S. 1987 in 300 Years of Gravitation, eds. S.W. Hawking and W. Israel, Cambridge 
University Press, Cambridge, 330. 



B. Allen: Cosmological sources 45 

Thorne, K. S. 1995 in Particle and Nuclear Astrophysics and Cosmology in the Next Millen- 
nium, eds. E. W. Kolb and R. Peccei, World Scientific, 160. 

Thorsett, S. E., Dewey, R. J. 1996 Phys. Rev. D 53, 3468. 

TSUBONO, K. 1995 in Gravitational Wave Experiments, proceedings of the Edoardo Amaldi 
Conference, World Scientific, 112. 

Turner, M. S. 1996, Detectability of inBation-produced gravitational waves, preprint 



astro- 



pli/9607066, submitted to Physical Review Letters. 
ViLENKiN, A. 1985 Phys. Rep. 121, 263. 

ViLENKiN, A., Shellard, S. 1994 Cosmic Strings and other Topological Defects, Cambridge 
University Press, Cambridge. 

Weinberg, S. 1972 Gravitation and Cosmology: Principles and Applications of the General 
Theory of Relativity, John Wiley & Sons, New York, 513. 

Wright, E. L., et al. 1992 Astrophys. J. 396, L13. 



